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ABSTRACT 


Maniscalco, James Andrew. Ph.D., Purdue University, 

May 1973. AN ANALYTICAL TREATMENT OF PHOTON TRANSPORT IN 
MumerAYERED MEDIA. Major Prcfessor: H. E. Hungerford. 

In this thesis the transport of photons in two-layered 
Slabs is examined, and then evaluated by analytically ex- 
pressing the two-medium buildup factor in terms of single 
medium buildup factors and parameters which in turn can be 
characterized as functions of atomic number. The word 
analytical is used to indicate that the parameters in this 
formulation are not empirical but rather they are directly 
related to the basic physics of transport phenomena in two- 
layered slabs. 

In applying the two-medium buildup factor formulation, 
we assume that the single medium buildup factors and albedos 
are available and then seek to determine the two-medium 
parameters. These quantities have been evaluated by using a 
one-dimensional discrete ordinates transport code (ANISN) to 
calculate the attenuation of photons in two-layered slabs 
of water, iron, and lead. Independent verification of some 
of the ANISN results was made using the stochastic transport 
code PUGTII. The materials water, iron, and lead were chosen 


because they are commonly used together in pair combinations 





XXV 


and because they represent the light, intermediate, and 
heavy elements respectively. 

In the process of tabulating and plotting the results, 
it was observed that the two-medium parameters are slowly 
varying functions of the first medium and smoothly varying 
functions of the second medium; moreover, a detailed examin- 
ation of the plotted results (Appendix D) reveals that the 
curves for each of the two-medium parameters have the same 
shape {i.e., functional dependence) for different first- 
layer thicknesses and/or materials. Because of this func- 
tional behavior, these two-medium parameters have been 
accurately extrapolated to a large number of two-material 
combinations by simple Minest interpolation on atomic num- 
ber. 

Energy and dose buildup factors for more than 15 differ- 
ent two-material combinations have been calculated with the 
two-medium buildup factor formulation and compared to cor- 
responding ANISN results. These comparisons clearly indi- 
cate that the errors resulting from obtaining buildup factors 


meen this formulation are less than 5%. 











T. INTRODUCTION 


1.1. Statement of the Problem 

The advent of the nuclear age has resulted in the de- 
velopment of sophisticated nuclear devices and systems in 
two vital areas of our material well-being: (1) national 
security; and (2) energy. The weapons and nuclear power 
reactors conceived and developed in the many military and 
civilian programs over the past quarter of a century have 
been possible only because of our advancement in the under- 
Standing of nuclear radiation transport processes within ma- 
terial media. As part of this theoretical undertaking, a 
considerable amount of effort has been devoted to developing 
methods of calculating the deep penetrations of photons in 
material media. These radiation transport calculations 
usually involve the solution of the Boltzmann transport equa- 
mon in One of its various forms. 

miessboltzmann transport equation is a statistical- 
mechanical statement of the balance and flow of radiation in 
terms of position, energy, direction, and time. As such, it 
is an integro-differential equation which can not be solved 
Meeaecly fOr any but a few highly restricted cases; however, 
numerous methods have neon developed and used to provide 


approximate solutions. They include techniques utilizing: 





Ime sile@cessive imtegrations over trial solutions, 
2. Polynomial expansions such as spherical harmonics, 
mee Fourier transforms, 


pee 6©Oitffusion theory, 


5. Variational techniques and perturbation theory, 
6. Expansions in terms of moments, 
fee invariant imbedding, 


8. Discrete ordinates S, theory, and 


N 

Memeo eocihwastic Solutions (Lattice Transport and Monte 

Barlo). 
Some of these methods most applicable to photon transport 
(the topic of interest in this investigation) will be described 
in the next chapter. 

It has become customary to present the results of photon 
transport theory calculations in the form of a buildup fac- 
tor. The buildup factor is defined as the ratio of the de- 
Sired quantity which is a characteristic of the total gamma- 
ray flux density to a similar quantity which is a character- 
istic of the unscattered flux density. For example, the 
energy buildup factor is defined as the ratio of the total 
energy flux (direct plus scattered) to the unscattered energy 
flux (direct). The principal advantage of expressing the 
results of photon attenuation calculations in the form of 
buildup factors is that these quantities are rather slowly 
varying Fue cei of attenuation distance, photon source 


energy, and average atomic composition. 





The penetration of gamma-rays in an infinite, homogeneous 
medium has been explored in detail using the moments meth- 
mel 23.4] a technique particularly suited to photon trans- 
fort in infinite media. Goldstein and alee. have used 
this method to compile an extensive collection of buildup 
factor data which for 20 years has been and is still con- 
mmered to be the primary source of buildup factor data for 
meee 1SOtropic and plane monodirectional sources in infinite 
fear. Of more practical interest to shielders is the trans- 
mission of photons through finite slabs of material. Much 
work has been done on single (one-layer) slabs, 12°0?/ 28] in- 
Searing calculations of reflection coefficients and transmis- 
$i0n buildup factors for several different source distribu- 
mes (isotropic, perpendicular, and oblique incidence). 

A small amount of work has been done on two-layer 
pet? 10511412) However the object of such studies has 
been to derive corrections to single medium buildup factors 


and/or empirical Formulast!3514.15.16] 


which work only in a 
given set of circumstances. Wo investigation that the author 
1S aware of has focused on expressing the two-medium buildup 
factor in terms of general parameters which are directly re- 
lated to the basic physics of transport phenomena in such 
layers. 

The objective of this thesis 1s to investigate tne trans- 
port of gamma-rays in two-layer slabs and to develop an 


analytical expression for two-medium buildup factors in terms 


of single-medium buildup factors and parameters which, in 





mite Can be Characterized as a function of atomic number. 
The reason for expressing the two-layer results in terms of 
one-layer results is that one-layer results are readily 
available in the open literature, and two-layer buildup 
factors so expressed will be of great practical benefit. 

The present theoretical investigation of photon transport 

in two-layered slabs was conducted by first thoroughly re- 
viewing the literature for all existent one- and two-layer 
meeaup and transmission information. As a result of this 
extensive search, it becane evident that the problem of two- 
layer buildup factors could best be handled by utilizing 

Seema iscrete ordinates and stochastic transport techniques 
to calculate, as accurately as possible, the attenuation of 
gamma-rays in two-layered material combinations of water, 
7ron, and lead. These elements were chosen for the study 
because they are commonly used together in pair combinations, 
and they represent the light, intermediate, and heavy elements 
respectively. The results of these calculations were then 
used to determine a set of parameters required in the formu- 
lation of the two-material buildup factors. 

As has already been pointed out, a number of computa- 
tional techniques are available for determining the transport 
Of gamma-rays in two-layered media. However, these tech- 
niques usually involve the extensive use of a large computing 
facility, and when a large number of calculations on complex 
Shielding configurations are to be undertaken, this method 


1s unduly slow, awkward, and costly. We looked, therefore, 





fereways to Simplify the approach for others, performing 

in advance (as it were) the long computer calculations for 
them. The unique feature of the analytical formulation 
derived in this work is that two-material buildup factors 

can be easily hand-calculated, using fairly simple formulas 
for most of the two-material combinations of interest to 
Smmera designers. In fact, the major limitation on the 
number of two-material combinations which can be investigated 


1s the lack of one-medium data. 





CE Hes RANSE ORT OF GAMMA RAYS: 


REVIEW OF THE LITERATURE AND METHODS 


The attenuation of gamma-rays in material media 1s rela- 
tively simple to predict when each interaction process 
meoerrsion) results in the disappearance of a photon, i.e., 
the medium is purely absorptive. For example, consider the 
case of a monodirectional beam of gamma-rays of energy E In 
Such an absorptive medium; then, the change in the flux in- 


memo ty of the beam dI(x), as it traverses the medium, is 


given by 
d] MeV 
ped = y 1(x) eC, ai 
cm~-sec 
where the constant of proportionality y is known as the 


linear attenuation coefficient. The solution of equation 
(2.1) is the well-known Lambert's law of absorption, 


Ti = a es (2.2) 


2 
cm -sec 
Unfortunately, Lambert's Law is not an accurate description 
Of photon attenuation in material media. First, the law 
assumes that the medium is purely absorptive and does not 


account for Scattering. Second, the law does not account for 





the energy change which results from photons of lower energy 
mentributing to the flux at x. 

Lambert's law can be corrected to account for scattered 
megtacion by using a multiplicative factor called the 
mmercup factor. the term “buildup factor’ was first intro- 
duced during the Manhattan Project. Its usefulness stems 
from the fact that radiation quantities resulting from un- 
Scattered photons, i.e., photons which eimerge without suf - 
fering any collisions, are relatively simple to calculate. 
There are many different quantities of interest in shield 
G@esign (particle flux, energy flux, dose rate, etc.) and 
Memean define a buildup factor for each of them. To show 
how this can be done, let the superscripts u and s refer to 
unscattered and scattered photons respectively and let the 
Guantity of interest be represented by the functional F{I); 


then, the buildup factor for the quantity F(I) is defined 
i 


aelest ch (een) 


In words, equation (2.3) defines the buildup factor as the 
mero Of a desired quantity which is a characteristic of the 
total flux density (unscattered plus scattered) to a similar 
met ity which is a characteristic of the unscattered flux 
Mmemoity. From this definition, it is obvious that the build- 
up factor concept is meaningless when the unscattered quan- 


Mmecy does not exist. 





Buildup factors of interest in shield design include: 


1. The energy buildup factor B, (with F(I) represent- 


E 
ing the integrated energy flux) defined as 


Ie ie} de 
Be (x) ” 7 7 a ; (245) 
iF: re {(x,£) de 
The total energy flux I(x,E&) is the sum of the scattered 


and unscattered contributions. That is, 


U S 


Meee La e(x,e) + 17 (x,E) (225) 


Hence, 


Be (x) =) + ———__—_—__.. cGy 


yee yhe dose buildup factor By With F(I) representing 


the dose-rate in air) defined as 


S 
es Geet (x.E) dé | a 


eee) 1 (x46) de 


By (x) 


where a (E) is the energy absorption coefficient of air. 


3. The energy absorption buildup factor By, wath 11) 
representing the energy density absorbed by the medium 


from the photons in unit time) defined as 


fare eel  (x,E) de 


Rex | + 
A Se ie (E) ite ieee) de 





where Or (E) is the energy absorption coefficient for the 
medium. It should be mentioned that in all of the defini- 
tions cited, the functional F depends only on the energy 
flux; therefore, for a monoenergetic source, Fei) is given 
mea Constant, Fins times the integrated, unscattered energy 
flux ex). In these cases, the buildup factor is simply 


defined by 


F(I) = BL F. I(x) . (229) 


Aoeunansport- Theory 


The fundamental equation governing the transport of 
gamma rays in material media is an integro-differential 
equation Known as the Boltzmann transport equation. ia hS 
basically a bookkeeping statement which accounts for photon 
increases and decreases in a given increment of space, 
angle, energy, and time. For shield analyses, a steady 
State is normally assumed and the time variable is not re- 
quired. Briefly, the equation of continuity is applied to a 
S1x-dimensional phase space consisting of three spatial 
coordinates r, two directional coordinates {) and one energy 
coordinate E; and the resulting integro-differential form 


of the Boltzmann transport equation is given as 


> 


We 6 (r.£.9) + vo(r.£,%) = S(r,E,%) 


+ om ee oir, —’, 8') de‘ do’ 
Ee) © 
2. 10) 
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In the above equation, the quantity »w is the linear attenu- 
ation coefficient and the angular flux o(r,£, 9) is defined 


so that 


o(r,£,8) dE dd (2.198 


gives the number of photons at r, with energy E in the 
range dE, going in the direction QO within the element of 
solid angle do, which cross in unit time a unit differential 
element of area whose normal is in the direction %. The 
Various terms in equation (2.10) are defined as follows: 
V- 8 d is the loss rate due to leakage from the phase 
space cell (dr, dE, d%), 

uo is the loss rate due to collisions in the phase space 
cell which either absorb the photons or change their eneray 
and direction, 

S 1s the independent generating source of photons, and 


GEE oo) o(r, E'S) dE* de (222) 


S 


> 
Bo) 
represents the rate at which photons are scattered into the 
phase space cell. The transport kernel, ie Eo 0), 
represents the differential macroscopic cross section for 
Scattering from the direction Q' to 2 and from the energy 
F* to — per unit solid angle and energy range. For Compton 


Scattered photons, the transport kernel is given by the 


well-known Klein-Nishina formula. This formula takes a 





particularly simple form if the energy dependence is ex- 


pressed in terms of photon wavelength in Compton units, 
AO z 51] 

A{( Compton Waits ) = E(MevV) = E(MeV) : (2.133 
In these units, the transport kernel is 

Mr '>),0'sa)=5 r : aty® = oe Ze a) 

a : ; 2 e N ae d 

meee tee 4° >i - 3-8) , (2.14) 

where n is the electron density and Yo is the classical 
meus Of an electron (2.817 x 10°! Sli rne Dirac delta 


function 6 accounts for the fact that wavelength (energy) 


and angle changes are related by the Compton shift equation 
MAho- AeT- te Gg. eral) 


In shielding problems where absorbed energy density 
and biological! dose are the quantities of interest, it is 
convenient to work with the angular energy flux which is 


defined as 


1(r,E,%) = E - o(r,£,@). (2.16) 
* One unit of Compton wavelength = h/m.C = .02425°A, 


where h is Planck's constant, m_ is the rest mass of an 
electron, and C is the velocity of light. 





Wee 


This preference arises from the fact that absorbed energy 
density and biological dose are tiore nearly proportional 
to energy flux than particle flux. The transport equation 
is written in terms of I(#,4,%) by multiplying equation 


meel0) by E and noting that 


3y)/) 


dE = oe ob a (2s 
A 
fe Sirek.o ee? ae (2ae)) 
1(r,—£,8) = 1(r,n,9) (2.19) 
and 
of: > ,¢ +, > 
EO = a 2 (6 ASA, 820) (2.20) 
then, 
mee T(r,4,9) + vw I(r,d,0) = S(r,2,98) 
q |, am eee ir 9. 0') di’ dp'.(2.21) 


A direct, analytical solution of the Boltzmann trans- 
port equation is feasible only for a very small number of 
Simple and highly idealistic problems. There are, however, 
many methods available which provide approximate solutions 
to the transport equation, several of which are described 
below. In general, methods utilizing numerical techniques 
with high speed computers provide more accurate solutions 
than methods which employ simplifying assumptions to alter 


the transport equation such that an analytical solution is 





obtainable. An example of the latter type is the elementary 
one-speed diffusion equation of reactor physics which can 
be derived from the transport equation by invoking the 
following assumptions: 
1. Scattering is isotropic in the laboratory co- 
ordinate system. 
2. There are no strong absorbers in the medium 
nes” De eee: 
3. The region of interest is two or more mean free 
paths from any intense sources or boundaries. 
The existence of these limiting assumptions (especially the 
first assumption) is a clear indication that diffusion 
theory is inadequate for problems involving photon trans- 
port because the scattering of photons is not isotropic at 


any energy. 


2.1.1. The Method of Successive Scattering 
A relatively simple and straightforward approach for 
calculating the contribution of scattered photons to the 
penetration of gamma-rays is found in the method of Ssuc- 


Oe] 


Meosive scatterings. This method is a@ semianalytica! 
technique which considers the total photon flux to be the 
sum of the unscattered flux, the once scattered flux, the 


twice scattered flux, etc. In notational form 


eee) 61x) (2.22) 





where >. (x) ig the flux at x which results from the ith 
scattering. As previously mentioned calculations for the 
unscattered flux are relatively simple to perform; however, 
this is not so for the scattered flux. In the successive 
Scattering method, the unscattered flux provides the 
collision density for first scattering which is used as a 
SOurce in calculating the once scattered flux. The process 
is then repeated with an integral recursion formula which 
gives a relationship between 247 'x) and >. (x). In other 
words, the higher-order scattered fluxes are expressed as 
integral operators of the next lower-order scattered flux. 
After the first scattering, the mathematics becomes quite 
involved and there is a tendency to limit this type of 
calculation to only one or two orders of ecg. |. 
However, this simplistic approach yields unreliable results 
for deep penetration problems where many scatterings con- 
tribute to the flux and the once and twice scattered con- 
tributions are only a small part of the total flux. 

The method of successive scattering was first used by 


[18,19] to calculate the transmission 


Peebles and Plesset 
Seemonodirectional gamma-rays through finite slabs of iron 
and lead. The error in these calculations is estimated to 
be as high as 20% for slab thicknesses up to 20 mean free 
paths. A variation of the method of successive scattering 


[20] 


has been developed by Peebles and elaborated upon Dy 


temush.-2!! ieweonststs of considering the photon 





transmission through a thick slab to be composed of a suc- 
cession of transmissions through thin slabs which are taken 
thin enough to make the contributions cf photons scattered 
three or more times negligible. Gamma-ray transmission 
in the thick slab is then obtained by successive matrix 


multiplication of the thin slab results. 


elec sine Spnerical Harmonics Method. 

The method of spherical harmonics, as applied to the 
transport equation, consists of expressing the various an- 
gular dependent terms in a complete set of elementary 
functions such as a series of polynomials. In three- 
dimensional geometry, spherical harmonics are the logical 
choice; however, these reduce to Legendre polynomials for 
one-dimensional geometry with azimuthal symmetry. When 
the polynomial expansions are substituted into the trans- 
port equation, it separates into an infinite set of coupled 
differential equations which no longer involve the direc- 
tional variables. Practical methods of solution require 
that the infinite set be limited to a finite number of 
coupled equations, and this is achieved by truncating the 


polynomial expansions to some finite number of terms. For 


“Although not widely used today, this method embodies 
concepts which are used in other techniques aescribed in 
this Chapter. It has therefore been included for peda- 
gogical reasons rather than its potential for problem sol- 
Pang. 
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example, a i (or nth order) approximation represents a 
solution in which the polynomial expansions are limited to 
n+] terms. 

In genera], the degree of anistropic scattering dic- 
tates the order of expansion required to yield accurate 
results. The Py approximation is equivalent to diffusion 
theory which, as previously indicated, provides adequate 
results only when the scattered flux is nearly isotropic. 
It has been shown that a P approximation provides a vast 
improvement over diffusion aay ed and that the ‘ost 
severe anisotropy resulting from Compton scattered photons 
Meee nergies Jess than 10 MeV can be adequately resolved 
with a is eee cmeriens tos “atin ag bee has used the 
Spherical harmonics method to investigate gamma-ray heating 
in multi-slab geometry. He successfully calculated the 
Spatial distribution of the energy flux spectra out to ten 


mean free paths with a P approximation. 


2.1.3. The Moments Method 
eo] 


The moments-method of Spencer and Fano is a semi- 
analytical technique which solves the Boltzmann transport 
equation in infinite geometry. With this method, as in 

the spherical harmonics method, the angular flux is expand- 
ed in terms of Legendre polynomials. “ In addition, spatial 
moments of the Legendre coefficients, Li xsE), which serve 


to separate energy and spatial dependence, are introduced 


by the definitions 





ie 


x,E) x” dx es 


3 





Eee I 


m! a! 


= OO 


where Uy is the linear attenuation coefficient at the source 
energy EO: In terms of these moments, the transport equa- 
tion reduces to a set of interlinked integral equations of 


the Volterra cpp eld 


with energy as the only independent 
variable. The spatial moments are straightforwardly eval- 
uated by numerically integrating these integral equations. 
meeeonould be emphasized that the calculations to this point 
are rigorous because no assumptions have been made. The 
major source of error in this method results from recon- 
Structing the spatial dependence of the flux with only a 
finite number of known moments. Usually six to eight 
properly chosen moments are sufficient to reconstruct the 
spatial flux with reasonable accuracy.* 

The method of moments is generally considered to be 
the technique of choice for problems involving photon trans- 
port in infinite media. Goldstein and ieee have used 
the moments method to investigate in detail the penetration 
of monoenergetic gamma-rays in infinite homogeneous media, 
and they have used their results to develop an extensive 


compilation of energy flux spectra and buildup factors for 


“The interested reader is referred to Reference [1] for 
a detailed description of two methods used to reconstruct 
the spatial flux The polynimial expansion method and the 
method of undetermined parameters. 





infinite plane and point isotropic sources. Advantages of 


raya 


the moments method, as stated by Golstein, include: 

1. A detailed energy spectrum of the flux can be 
obtained. 

2. The required computing time for penetrations up 
to 20 mean free paths is relatively short. 

3. Most photon source configurations of interest 
are amenable to the method. (Spencer and 
cEve| n echoes have successfully employed the 
method to calculate the slant penetrations of 
photons in water and concrete.) 

4. Angular distributions of the flux may b2 obtained. 

[30] 


(Certaine used the moments method to calculate 
angular distributions from plane monoenergetic 
sources.) 

The most serious limitation of the moments method is 
that it 1s applicable only to infinite homogeneous media. 
Another disadvantage, less restrictive but still noteworthy, 
involves the difficulty in reconstructing and determining 
the flux behavior near the source (less than one mean free 
path). Because of this, the method of moments is not a 


good technique for evaluating reflection coefficients (al- 


bedos). 


evaerl> 4 has developed a method to alleviate the 


mereticuity for point isotropic sources. 
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2.1.4. The Method of Invariant Imbedding 

The invariant imbedding technique was originally devel- 
oped and used to calculate the reflection of diffuse light 
from a stellar atmosphere by the astrophysist Ambrazum- 
fan. 632331 Bellman, Kalaba, and Wm have extended 
the method into a general approach, and shown that it can 
be applied to a much broader class of problems which in- 
cludes neutron and gamma-ray transport. As apolied to radi- 
ation, the invariant embedding formulation is not another 
method for solving the transport equation but rather a new 
fundamental approach to the transpcrt cf radiation in ma- 
terial media. 

Unlike the Boltzmann equation, the dependent variables 
of the invariant embedding formulation are reflection and 
transmission functions. Each of these parameters is de- 
fined by and thus satisfies a non-linear integro-differential] 
equation which cannot be solved analytically. Hence, all] 
of the solutions are numerical in nature. The reflection 
and transmission equations represent problems of the 
"initial-value" type. In contrast, the Boltzmann equation, 
as applied to photon transport in slab geometries, forms 
“boundary-value" type problems. In general, initial-value 
problems are more amenable to numerical solution with a 
computer. As a consequence, the invariant embedding formu- 
lation has a computational advantage over the Boltzmann 


approach. 
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Mizuta and scien os have applied the method of 
Invariant imbedding to the problem of photon transport in 
Slabs and demonstrated that it is competitive in accuracy 
to other established methods such as Monte Carlo and the 
method of moments. Compared to stochastic methods, in- 
Variant imbedding has the advantage that it requires less 
menouting time. Compared to the moments method, it has 
the advantage of being applicable to slab geometries. The 
G@isadvantages of the method are that it is difficult to 
apply to other than slab geometry and that it is inefficient 
for thin shields. (Thin shields require almost as much com- 


puting time as thick shields.) 


Vee The Discrete Ordinates S,, Method 


N 
om C. Mieke! 2 in 1943 was the first to suggest the use 
of discrete ordinates (or discrete directions) in transport 
theory. In the years that followed, S. Siamese lie eee 
further developed the method and applied it to the study of 
radiation transport in stellar atmospheres. Early applica- 
tions of the discrete ordinates method were primarily 
limited to the transport of isotropically scattered radia- 


tion in slab geometries. alow! with 


Some success, 
applied discrete ordinates to curvilinear geometries by 


uSing an “Sy SourOoOximation. In the Sy version, the angular 


dependence of the flux is approximated by a series of con- 


nected straight line segments; hence, the letter S signifies 


Segments and the subscript N indicates the number of segments 
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used to construct the angular dependence. It was subse- 
quently recognized by Carlson and tagprean la! that the 

Sy approximation was only a special case of a more general 
formulation and this led to the develonment of the discrete 
ordinates Sy method. 

The discrete ordinates Sy method, as presently employed, 
iS a means of numerically solving the Boltzmann transport 
equation with finite difference equations. The finite dif- 
ference equations are formulated as flow balances for dif- 
ferential phase space cells (dV dE ds). The more recent 
versions of the method accurately treat anisotropic scat- 
tering by using a Legendre series expansion to describe the 
scattering cross section, og GE ey tee eo The fundamental 
assumption in the discrete ordinates method is that the 
inscattering integral can be approximated by a quadrature 
Scheme which evaluates the integral as the sum of a dis- 
crete distribution. In one-dimensional geometries a Gauss- 


bad 


Legendre quadrature formula may be used to provide the 
discrete distribution, with discrete direction cosines cor- 
responding to the Gaussian Zeroes. 

The discrete ordinates Sy method has been shown to be 
a versatile and powerful method for accurately solving the 
transport equation in both one-dimensional geometries 
(DTF ry, 642] antsnt43J) and two-dimensional geometries 
(pot e444). It has been successfully applied to neutron and 


aaa 47 } 


gamma-ray deep penetration calculations as well 
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(48 ] The results 


as to general core criticality problems. 
of these caiculations and many others clearly indicate that 
the discrete ordinates method, as applied to shielding prob- 
lems, has the following advantages: 

1. Depending somewhat on the sophistication desired, 


the S,, calculations are easy to prepare. 


N 

meee ine metnod 1s not stochastic and flux errors 
at deep penetration are systematic rather than 
Seateistical. 

3. Secondary gammas may be calculated by the same 
method, either as a second calculation or simul- 
taneously with the neutrons. The gamma yield 
distribution may be dependent on the neutron 
energy at capture. 

4. The neutron energy range from highest fission ener- 
gies to thermal, including upscatter, may be cal- 
culated by the same method. 

5. The one-dimensional calculations are much faster 
(in computer time) than similar Monte Carlo cal- 
culations. In two dimensions the type of problem 


and the desired answers determine whether Sy or 


Monte Carlo is better. 


“These saan ta ges were compiled from a listing given by 
meer. Mynatt : 
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The major disadvantage of the method is that it is not 
applicable to three-dimensional geometries. 

The discrete ordinates method is one of the methods 
used in this work to calculate the penetration of photons 
in two-layered slabs; therefore, a derivation is given in 
Appendix A which establishes the equivalence of the one- 
dimensional discrete ordinates equations with the transport 
equation and which also shows how the Gaussian quadrature is 


employed to evaluate the inscattering integral. 


Zales otocnastic Methods 

Many physical processes including the emission and 
transport of radiation are probabilistic in nature. As a 
result, many types of problems in mathematics and physics 
can be solved by random sampling or stochastic techniques. 
The Monte Carlo method, as applied to radiation transport, 
represents a system of probabilistic solutions in which 
probabilities relating to microscopic emission and interaction 
processes are used to predict macroscopic phenomena such as 
flux and dose-rate. Simply stated, the Monte Carlo approach 
consists of constructing a series of analogue experiments 
(histories) which mathematically simulates the emission and 
transport of particles, and then analyzing the outcomes of 
a large number of histories to obtain relevant quantities. 
Each history (or experinent) includes the generation of a 
particle, its random walk through the medium, and its death 


which results when the particle is absorbed, exceeds some 
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physical energy or time boundary, or when its importance 
becomes insignificant. 

The method of random sampling was first applied to the 
Solution of radiation transport problems by Ulam and Von 


[oR imiecmma-io+0 Ss in the early stages of 


Neumann 
development, the Monte Carlo method was based on direct 
Simulation of the transport problem and was therefore found 
to be prohibitively castly and time consuming in calculat- 
ing deep penetrations with reasonable accuracy. In the 
years that followed, several techniques were developed to 
reduce the variance associated with the statistical behavior 
of a particle history. Some of these variance reducing (or 


[52,53] include: 


biasing) techniques, as described by Kahn, 
1. Importance Sampling: The selection of events 

(source parameters, path lengths, etc.) or omis- 
sion of events (absorption) is based on an altered 
rrobability distribution which allows the more 
important events to be sampled more frequently. 
The exponential transform is one of the more widely 
used means of importance sampling because it is 
easier to alter selection of path length. Particle 
weighting is another widely used importance samp- 
ling technique in which absorption is not allowed 


but accounted for by reducing the weight of a 


Scattered particle. 
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2. Use of Expected Value: A portion of the random 
walk process is replaced by its average value. 

3. Russian Roulette and Splitting: Unimportant his- 
tories (i.e., those which are unlikely to contribute 
to the answer) are terminated with some probability 
while important histories are multiplied to give 
several semi-independent estimates. 

In general, Monte Carlo is not the technique of choice 
for one-dimensional steady-state transport problems because 
other methods such as discrete ordinates are usually much 
faster (in computer time) and just as accurate. In two- 
dimensions the type of problem (deep penetration, irregular 
geometry, ducts) eeetemn Ties whether Monte Carlo or discrete 
ordinates is better; but, there 1s no competitor to Monte 
Carlo for obtaining rigorous solutions to three-dimensional 
Or time-dependent transport problems. 

miewbattice Model Concept of Stochastic Transport was 


developed and successfully applied to neutron transport by 


q, 154] It differs from conventional Monte 


Meee. Hungerfor 
Carlo methods in that the random walks of particles are 
constructed by sampling their motion at predetermined points 
uSing a direct simulation of the transport process. In 

the “Lattice Model Concept" space can be envisioned as being 
made up of a large number of small cubes whose sides are 


equal to a distance called the unit lattice distance. The 


method of sampling particle motion within this framework is 
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(541 


Meet described by H. E. Hungerfor in the following 


excerpt: 


If we now restrict the locations of the sampling 
points to be only at the corners of the cubes, and 
nowhere else, we have already chosen our allowed 
directions of emission, namely, these directions 
represented by lines which can be drawn from any 
corner of the unit cube to any other, neglecting 
the repeats. There are altogether 26 allowed 
directions of motion, outward from any given source 
point (considered at the origin)--6 along the 6 
axial directions, 12 in-plane diagonal directions 
(the planes being defined by any 2 of 3 axes), and 
8 skew diagonal directions (to opposite corners of 
the cube). Actually, to show all these directions 
emanating from one point, one has to make a larger 
cube from 8 unit cubes and place the source at the 
central point common to all cubes. These 8 cubes 
correspond to the 8 octants of a sphere. The larger 
cube [may be] referred to as the representative cube. 


[55 ] 


A. Razani applied the lattice concept to gamma-ray 
transport by developing a stochastic gamma-ray transport 
code called PUGT I. Another version of the code, PUGT II, 
which uses particle weighting will be used in this work to 
calculate the penetration of photons in two-layered slabs, 
and thereby provide independent verification of the discrete 
Ordinates results; therefore, a derivation which relates 


the transport of photons in the lattice model to the inte- 


gral form of the transport equation is given in Appendix B. 
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3. THEORY OF PHOTON TRANSMISSION AND 
BUILDUP IN A TWO-LAYER SLAB 


3.1. Analytical Formulation of the 
Two-Medium Buildup Factor 


Given two slabs of different materials and a plane, 
monoenergetic gamma-radiation source incident on one of them, 
the problem under consideration is to develop a relatively 
Simple expression which will predict the buildup and there- 
fore the transmission of radiation through the slabs, with- 
out having to resort to expensive and time-consuming com- 
puter calculations. The choice of the radiation quantity 
selected as the prime Variable is the angular energy flux 
density, Mibaae 0) in MeV per cmé per second per steradian 


and per unit energy, where 


Mapes =F = 6(b- £2) 


ikea | 


and (b,,E, 5) is the angular photon flux density. (The num- 
ber of photons of energy E at b, Teossmooetnrougn a unit of 
area in the direction 0), per unit time, per steradian, and 


per unit energy.) The analytical formulation to be 


“Quantities used in this development are defined in the 
Nomenclature preceding the text. 
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developed is based on the buildup factor concepts, as given 
in Chapter 2, and the geometrical configuration shown in 
Figure 3.1. More rigorously stated and referring to Figure 
3.1, the objective of the formulation to be developed is 

to determine the buildup factor which when used with Lam- 
bert's Law correctly predicts the energy flux density at 
the detector (point C) resulting from a plane, monoener- 
getic, monodirectional photon source at point A. The two- 
medium buildup factor, as derived, will be a function of 
Guantities which characterize the angular energy flux at 
the interface between the two slabs. AS nearly as possible 
the formulation for two-layer buildups will follow that for 
One layer. The Single-medium buildup factors previously 
defined contained two terms, one describing the scattered 
energy flux and, the other describing the unscattered ener- 
gy flux. The same scheme will be followed here. In addi- 
tion, a distinction will be made between photons which 
Scatter in the first layer and those that scatter in the 
second layer. In order that the reader fully understand 
the terminology and import of this development, a series 

of terms used herein is defined. The superscript s is used 
to refer to the scattered radiation and the superscript u 
refers to the unscattered radiation. These superscripts 
may be used alone or in combination. Thus: 


1. I° is the scattered energy flux, 
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Figure 3.1. Geometrical Configuration for a Two-Layer Slab. 
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2. I is the unscattered energy flux, 


fm 1°" 


refers to the energy flux which has scattered 
in the first layer but not in the second layer, 

4. °° refers to the energy flux which doesn't scat- 
ter in the first layer, but scatters in the second 
layer, and 

5. I>> refers to the energy flux which scatters in 
both layers. 

These and other terms will be more fully described as the 
detailed development unfolds. - 


The scattered energy flux at the detector (point C) 


may now be written as the sum of three components 


SS MeV 
) + ] (b, +b.) ? 5 
cm sec 


) + 14> (b +b 


aS 
= |  (b.+b )*b> 


S 
eb. +b +b, 


] 2) 


foe) 


where feb. +b.) is the scattered energy flux at the de- 
tector arising from the scattered energy flux at the inter- 
face which doesn't scatter again in traversing the second 
layer, ee by) is the scattered energy flux at the de- 
tector arising from the unscattered flux at the interface 
which scatters at least once in traversing the second layer, 
and me (b,+b,) is the scattered energy flux at the detector 
arising from the. scattered energy flux at the interface 
which scatters at least one more time in traversing the 


Second layer. Here, as in the rest of this work, the 
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quantity by denotes the number of mean free paths in the 


first slab at the source energy ES and b. 1s analogously 


defined for the second slab. 
b, = Uo (E) Xx, where My and Mo 
coefficients of the first and 


The relationship between 


That ) x and 


leS5 by 7 (E. : 
are the linear attenuation 
second media. 


the energy flux (Sade) 


and tne scattered angular energy flux at the interface, 


meee, E, <2), is given by 


f(b, +b 


>) = dE 


; An 


(3329 


In slab geometry the angular energy flux can be expressed 


aS 


S 
i ber cw) 
B(b,.£,8) = - mew ; eos 
cm -sec-MeV-STER. 
and the differential solid angle is given by 
dj) = w os Stendaians, (3.4) 
e 


where w = 


Stituting these relationships 


fh -(b +b 


cosé, and 0 and r are shown in Figure 3.1. 


1>(b,,E,w) e 


Sub- 


imeomeaquation (3.2) yields 


Oe ye dy MeV 
3 


fe CM =S ee 


f35) 





a2 


The integration variable may be changed to w by using the 


x5 i) Ree 


ferring again to Figure 3.1, one sees that w+Oas re oo, 


relations op = X5/W, r= X» tan6, and rdr = - 


and w > 1 as r>Q. Substitution of these variables and 


mts into equation (3.5) yields 


E 1 
-b 
ib, +b.) = dE ta E.¥) e dw 


0 0 


where b.(E) = No (E)x,. 

In proceeding further it will be necessary to expand 
the scattered, angular energy flux, 1° (b),E,m) in a series 
of Legendre polynomials. For a plane source, the Legendre 


Series expansion takes the form 


I> (b 





_ 22+] S 
; E,W) ~ ) Z 1) (b,>b,>E) Po (w) : (oe 
C= 


O 
It can be shown that beyond a certain number of terms ¢ = m, 
there 1s a negligible contribution to the scattered angular 
energy flux; therefore, the series is truncated at m where 

m indicates the order of the expansion. The terms Po(w) 
are the Legendre polynomials of order g, and the expansion 


coefficients are given by 


2 I5(b,,E,w) Po(w) dw Mev 


-] cm -sec-MeV 


I? 


g(b,>b,, (3.8) 


eit is anticipated Ehat a fifth order expansion will be 
Sufficient to describe I (b,,E,w) to very good accuracy. 
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The physical significance of the expansion coefficients can 


Memocertained by noting that for 2 = 0 and 2 = 1, these quan- 


tities represent the scattered energy flux spectrum and the 


Scattered energy current spectrum at the interface between the 


two layers. The higher g terms define moments of these. In 


addition, the scattered energy flux at the interface is relat- 


ed to the zeroeth order Legendre expansion coefficient by 


Fo 


S _ MeV 
I (b,.b,) = | I 
O 


2 
cm -sec 


= 


; E) dE 


(b,.b (3.9) 


tis 


Substituting the Legendre expansion for Goan EW ) into 


equation (3.6) gives 


O -b,(E)/w 


mM 
Mamipss>,) = y tl emTaChesh E) P. (Ww) e 


y tbo) 2 2 2 


2=0 O O 
MeV 


Cie Sec 


3m] 6D) 


For convenience in later analyzing and tabulating, it is de- 
Sirable to consider a two-medium parameter which is a slowly 
varying function of the thickness of the first slab. This 
quantity results when both sides of equation (3.10) are di- 
vided by the scattered energy flux at the interface as de- 


fined in equation (3.9). Thus, 





m E as 
O Gores Diag © ) -b,(E)/w 
eee) ee | cet | PS wy ee * 
n ie 2 2 is ee) Q 
Q=0 O O eaeZ 


et) 


dw 
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where 1°" (b +b.) is nothing more than eee o. | normalized 
to a scattered energy flux at the interface of ey cil. 
For convenience in carrying out the integration with respect 


to w, let b = bo(E). fete that b + b(E.). RT Soule t “us 


define a set of quantities 
] 


mM {b) = p(w) Peed w Gar > = 1,2,....m. (3.12) 


Substitution of M (db) inm-oreqdiiatiaon (3.11) yields 


m E S 
ole lope | cae aie 
Memiesb,) = jy Sgt eS yb) de. Gag 
n ie 2 2 eh b.) Q 
2=0 0 hier 


The evaluation of M,(b) is accomplished by comparing it to 


the exponential integrals which are defined as 


OO 
~V 
n- | e dv 
n n 
V 


(3. 14) 
The relationship between M,(b) and the exponential integrals 
1S Made clearer by substituting v = b/w into equation (3.12) 


So that 


(340169) 
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for ¢g = i, Pi) (b/v) Shey. anid 
M, (b) = E,(b), (31® 
for ¢ = 2, Po(b/v) = 4[3(b/v)* - 1], anc 
= . 
M.(b) = xL3 7a, = E,(b) J, ( Saes 
Eeec. 


Since kw) is a polynomial in w, M.(b) can be expressed 


Q 
in shorthand form as 


M(b) = PLEs,,(b)] - (3.19) 


Here the somewhat strange notation signifies that whenever 


w2 occurs in P tw), it is to be replaced by Es 42(b). For 
example, when g = 3, 
| 3 
P.(w) = 5L5 w - 3 wl , (32205) 
and 
ees _ 
M(b) = x15 E,(b) 3 E,(b) J foc ip 


With Mb) now expressed as an exponential integral 
function of energy, the final step in evaluating eet o bye 
as indicated in equation (3.13), involves integration with 
respect to energy. In performing this integration it is 
1° (b, bo, 
therefore, the energy flux ie(b 2b.) will be calculated in 


necessary to specify the expansion coefficients, Eps 


the next chapter after the expansion coefficients are de- 


termined. 
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The scattered energy flux at the detector which was 


previously expressed as 


MeV 
Gila sec 


= 1°4(b +b 


S 


us SS 
>) + [°°(b +b.) + | (b, +b.) 


] 


(3 eee) 


can be displayed in a more convenient form by defining a 
Scattered energy buildup factor, 2 (Soe, as the variable 
which gives the buildup of the scattered energy flux at the 


interface as it traverses the second medium. That is, 


SU SS 
I (b, +b.) ana (b,+b,) 
2 Su : 
I (b, +b.) 


ieee) 
+. ————____ = . ES 2389) 


Su 
I (b, +b.) 


With the scattered energy buildup factor defined in this 


manner, it follows that 


MeV 


cm =see 


Moe (bet b 


= > 
I°>(b,+b = £8 >) 


SU 
>) 7 ie (b,+b 


ie 2) i 2 


(3.24) 


Further simplification iS achieved by expressing "> in terms 
of the previously defined single-medium energy buildup fac- 


bor, that is 


) 
Mmibs+b,) = Ee | * [B.(b,)- 1] —pel , (3.28) 





oi 


where 
Tb, tb ) 
Be(b,) = 1 + Ee RAG) 
I (b, +b.) 
and 
= bet be) 
["(b,+b,) = E, e i oeeeetevi es (3.27) 


2 
Ell) “Ss Ele 


Substituting this expression for i into equation (3.24) 


yields 


(3.28) 


In keeping with the definitions for single-medium 
buildup factors, the two-medium energy buildup factor, 


Be (b,+b,), is defined as the ratio of total energy flux to 


unscattered energy flux. Specifically, in formula notation 


ee be p 


U 
ee baat >) 


tb | 2 


Be(b, +b,) Y, 


b, +b.) 


bei) 


Se (e273) 
“(ees oe) 


FE 6 e 
O 


Substituting the expression for ado by as given in 


memation (3.28), into equation (3.29) yieids 
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> S SU 
Paco (pb. ,b,) I- ({b,+b.) 
J Sb baa we n ee 


Ee 
0 


Equation (3.30) can be further simplified by noting that the 
energy buildup factor at the interface between the two layers 
iS given by 

S 

I (b),b,) 


Be (b),b,) = | + ates d (cies 


EY e 


where f(b, db.) is the scattered energy flux at the inter- 


face. mAus , 


S Su 
Be (b,,b,)[B,-(b,,b,)-1] (by +b.) 


B.(b, +b.) = Be(b.) - PY a 


This 1s the formulation which will be used to evaluate energy 
buildup factors for two-material combinations, other than 
those of water, iron, and lead which are calculated in the 
next chapter. [It will be shown in the discussion of results 
that the two-medium parameters BD D5) and (elon) are 
Slowly varying functions of b, and smoothly varying functions 
of bs which can be accurately evaluated for many two-material 
combinations by simple linear interpolation over atomic num- 
ber. 

The energy buildup factor at the interface between the 


two layers, B-(b, 4b Specdamturtener clarification so that 


a 
it can be evaluated from available single-medium parameters 


such as buildup factors and albedos. The scattered energy 
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mperat tne interface 1s composed of the following fluxes: 


f(b, ) - The scattered energy flux which results 
at the interface when the second medium 
1s replaced by a vacuum. 

i, ,b.) - The scattered energy flux at the inter- 


face which is reflected from the second 
layer. 


In notational form, 


S S MeV 
Sinise SEC 


Now, the scattered energy flux which is reflected can be 
expressed in terms of flux albedos as 

S a U S 

1.(b,3b>) = | (by) Ot 9) ¢(b,) OL (b),b,) <a » (3.34) 
where the unscattered energy flux albedo a 1s defined as 
the ratio of reflected energy flux to incident, unscattered 
energy flux; and the scattered energy flux albedo a 1s de- 
mined as the ratio of reflected energy flux to incident, 
scattered energy flux. Since the unscattered energy flux 

1s both monoenergetic and monodirectional, aU 


p 
to the energy flux albedos for normally incident sources 


is equivalent 


found in the open lq wapatenae Uae However, the scattered 
energy flux albedo is not dealt with in the open literature; 
therefore, it will be calculated and tabulated in this 


work. 
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Substituting equation (3.34) into equation (3.33) and 


dividing by the unscattered energy flux I (b,) yields 
= U > S 
B.(b,,b.) Be (b, ) + a4 (bo) + [B,(b, ) 1] ag(b,.b>), (32350) 


where Be (b, ) is the single-medium energy buildup factor for a 
normally incident source in slab geometry which has been 
evaluated in Reference [5]. It will be shown in the dis- 
cussion of rasults that the interface energy buildup fac- 

tor can be accurately evaluated for a large number of two- 


material combinations by using equation (3.35). 
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4. CALCULATIONS 


4.1. General 

The transport of monoenergetic gamma-rays normally 
incident on two-layered slabs of water, iron, and lead has 
been investigated by using both discrete ordinates and sto- 
chastic methods to solve the Boltzmann transport equation 
in slab geometry. Source energies of 1, 3, and 6 MeV were 
chosen so that the two-medium buildup factor formulation, 
as derived in the previous chapter, could later be verified 
Beeeene Primary photon interaction processes: The photoelec- 
tric effect, Compton scattering, and pair production. The 
production of annihilation radiation has also been included 
in these calculations, but other secondary interaction pro- 
cesses such as fluorescence radiation and bremsstrahlung 
have not. The consequences of not dealing with these photon 
interactions in the calculations will be considered and 
accounted for in the discussion of results. 

The discrete ordinates calculations were carried out 
with the Westinghouse version of a transport code called 
antsy, Loo] and the stochastic calculations were carried out 
with PUGT a gees a gamma-ray transport code developed at 
Purdue University. On the average, the PUGT II stochastic 


calculations required much more computer running time than 
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the ANISN discrete ordinates calculations, especially for 
slab configurations of more than three mean free paths. 
Therefore, the stochastic calculations were limited to two- 
layered slab configurations of up to five mean free paths 
and primarily used to verify the accuracy of the discrete 
Ordinates results. One of the major goals of this work 

will be to establish definitive accuracy limits for the two- 


layered buildup factor results. 


ome eCechaseirca Calculations 

Over 150 independent stochastic calculations have been 
carried out with the transport code PUGT II to evaluate the 
transmission and buildup of monoenergetic photons normally 
incident on two-layered slabs of water, iron, and lead. 
PUGT II (Purdue University Gamma-Ray Transport) is a sto- 
chastic computer code which is based on the "Lattice Model 
Concept’ as discussed in Chapter 2 and Appendix B. In brief, 
the code generates source gamma-rays and runs analogue ex- 
periments in which the photons are followed through a ma- 
terial medium until they exceed some physical, energy, or 
time boundary. The record of each photon from birth to 
Geath is called its history, and quantities of interest such 
as fluxes, currents, and albedos are obtained by summing 
the results of many thousands of histories. Reference (55) 


Should be consulted for a more detailed description of PUGT [J]. 





“A typical ANISN run for slab configurations up to 13 
mean free paths required approximately 250 seconds of time on 
the CDC 6500 computer at Purdue University while a PUGT II 
run for a slab configuration of 4 mean free paths required 
1300 seconds of computer time. 
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The results of these stochastic calculations for two- 
layered slabs of water, iron, and lead will be primarily 
used to determine the accuracy of the discrete ordinates 
calculations; therefore, it 1S imperative to estimate and 
limit the inherent statistical inaccuracies. For statistical 
purposes, the transmission of photons through a slab can be 
considered to be a series of independent Bernoulli trials 
where M, the fractional number of photons transmitted through 
the slab in N histories (trials), represents the probability 
that a photon will be transmitted through the slab on any 
given history. The standard deviation for this statistical 
quantity is then estimated by 


1/2 
eS PAE ay ea 1) 


N 


and the fractional deviation which indicates the percentage 
of deviation from the mean value of the fractional number 


of photons transmitted is given as 


F = 100 ye i .2) 


zn 


It is believed that a fractional deviation limit of 2% will 
provide a sound basis for determining the accuracy of the 
discrete ordinates calculations; therefore, the number of 
Nnistories for each of the stochastic problems in this work 
Was chosen large enough to limit the fractional deviation to 
within 2%. On the average 10,000, 20,000, and 30,000 


“The fractional deviation was less than 1.5% for problems 
1n which the second material was not lead. 
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histories were required for slab configurations totaling 2, 
3, and 4 mean free paths respectively. 

It should be noted that the statistical model, as 
presented above, does not account for the effects of biasing. 
Particle weighting is the only biasing technique used in the 
PUGT II code. This weighting procedure reduces the sta- 
tistical variance by analytically accounting for photoelectric 
absorption. However, the resulting decrease in the frac- 


tional deviation cannot be quantitatively determined. 


4.3. Discrete Ordinates Calculations 

Over 900 independent discrete ordinates calculations 
were carried out to calculate,as accurately as possible, 
the buildup and transmission of monoenergetic gamma-rays 
normally incident on two-layered slabs. The bulk of these 
calculations dealt with two-layered slabs of water, iron 
and lead, source energies of 1, 3, and 6 MeV, and shield 
thicknesses of up to 5 mean free paths of the first layer 
followed by up to 8 mean free paths of the second layer. 
The materials water, iron, and lead. were chosen because they 
are commonly used together in pair combinations, and they 
represent the light, cme cniaicer and heavy elements re- 
Spectively. The source energies of 1, 3, and 6 MeV were 
chosen in order to account for buildup and transmission char- 
acteristics resulting from all of the primary gamma-ray 
interaction processes: The photoelectric effect, pair pro- 


duction, and Compton scattering. 
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The discrete ordinates calculations in this work were 
carried out with the one-dimensional, discrete-ordinates 
code anrsn. £964 This computer program 1s written in 
FORTRAN IV and solves the one-dimensional Boltzmann transport 
equation in slab, cylindrical or spherical geometries. 

Theoretical considerations, as discussed in Appendix A, 
indicate that solution accuracy is highly dependent on the 
following parameters: 

mee Order of scattering oie 

2. Order of angular quadrature (Sy). 

meee Mesh interval spacing ( r). 

4. Cross section set used. 
The order of scattering ee represents the number of terms 
in the truncated Legendre series used to approximate the 
inscattering function. In general, the degree of anisotropic 
Scattering dictates the order of scattering required to give 
Mectrate results. A (P,) expansion was selected for the 
ANISN calculations in this work. 

The quadrature data sets used in the ANISN code con- 
Sist of direction cosines and associated weights which spec- 
My the discrete directions over which the inscattering 
function is evaluated. Symmetric sets corresponding to the 
Gaussian zeroes are given for various orders of angular 


[57] 


Quadrature .in the ANISN User's Manual. In general, for 
marsOtropic scattering, the order of angular quadrature 


(Sy) Should be at least twice as large as the order of 
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* 


meattering (P,). An (S,_) quadrature set in slab geometry 


QR 
was used for all of the discrete ordinates calculations in 
this work. 

The discrete ordinates derivation given in Appendix A 
emphasizes that solution accuracy requires adequate mesh 
Spacing. Stringent requirements to insure adequate mesh 
Spacing are presented in the ANISN User's onus ae The 
calculations in this thesis comply with these criteria in 
all cases. Specifically, the mesh interval spacings for 


these calculations were kept smaller than (4 i 


Eg 1s the largest total group cross section. 


, where 


ANISN fixed source calculations require the following 
group-averaged cross sections 
On - the photo electric and pair production ab- 
Sorption cross sections 


Oy [Uiecmitatalsattenwation cross section (total 
= absorption + scattering) 


Gq+g the ingroup scattering cross section 
Gees We= 15c,...,9-! - The Po down scattering cross 


sections. 
These group-averaged cross sections were generated with the 
ES.) ) 


Westinghouse version of a computer program called GAMLEG-W.: 


Kx 
With the absorption cross-section data and the group energy 


*k 
The quadrature requirements are treated in Appendix A. 


** 
The photoelectric and pair-production cross section 
data were qhtained from the tabulations of Storm and 
Israel. 
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bounds as input, this code performs the numerical integra- 
tions necessary to produce group-averaged absorption, scat- 
tering and total cross sections. fhe down scattering cross 
sections are determined by averaging the Klein-Nishina dif- 
ferential scattering cross section over the appropriate 
initial and final photon energy intervals. This version of 
GAMLEG does not account for the production of the annihila- 
tion radiation that results when positions created during 
pair-production events are destroyed. The code was there- 
fore modified to include the production of annhiliation 
radiation by treating this generation process as a pseudo- 
scattering event. In this modification, the ‘two (0.511MeV) 
annhiliation photons are assumed to be emitted isotropically 
at the same time and place that the pair-production event 
occurs, since the lifetime and range of a position are 
relatively very short in comparison to photon transport. 
With this assumption, the zeroth order down scattering cross 


Sections are modified according to 


; 


1+g 1 oeq = 1>g ° = 
 , (modified) oe ad ip 2 S(g p ) eas 3 ) 
where op is the average pair production cross section for 


group i, p is the index number of the energy group corres- 
ponding to 0.511MeV, and & represents the Dirac delta func- 
eon. 

Flux weighted 20-group cross section sets were used 


mor all of the Pe - 346 discrete ordinates calculations. 
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The 20-group energy structures used for the 1, 3, and 6MeV 
sources are tabulated in Tables C.59 - C.6]. The weighting 
fluxes were obtained from the infinite-medium energy flux 
Spectrum tabulations of Goldstein and Wilkins’! When the 
energy flux spectra, as calculated in this work, differed 
Significantly from the corresponding infinite medium tabu- 
lations, the cross sections were reweighted using the energy 
flux spectra obtained in the preliminary solutions. 

The following quantities were obtained directly from 


the output of the ANISN code: 


ib, +b,,£) - The scattered energy flux spectrum 
at the outer edge of the second slab. 

Be (b,+b,) - The two-medium energy buildup factor. 

rib, 3b.) - The Legendre expansion coefficients 


[mer ¢-1,2,...5 of the scattered angular energy flux 
at the interface, ees 
bo) The scattered energy flux at the 


interface. Note that, 


: 
O 


S = S 
I (b,,b,) = 15 (b,.b,,E) dE 
The two-medium energy buildup factors, BL (b,+b,), and the 
normalized scattered energy flux spectra, 
Map, +b,,£)/1°(b,+b,), are tabulated in Appendix C for three 


SOurce energies, six two-material configurations, and 


forty thickness combinations. The Legendre expansion 
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coefficients, Mere ESE), are used in the next section to 


71 
calculate the energy flux which scatters in the first layer 
Meo in the second layer while the scattered energy flux 
at the interface is used to determine the interface energy 

Syrdadup factor, B which in turn 1s used to evaluate 


p(b, »b,), 


the scattered energy flux albedo, a°(b,,b 


d line? >): 
4.4. Calculation of the Parameters in the 


Two-Medium Buildup Factor Formulation 


The two-medium buildup factor formulation which was 


derived in Chapter 3 is rewritten here as 


Beacper ie betb-4b,)-1]1."(b,+b.) 
Be(b,+b,) = 8-(b,) + a oe 
where 
dy .d) = B,(b,) + ay (bo) be a, (by »bo) ese) 


For organizing purposes, the quantities in this formulation 


are separated into two categories: 


1. Single-medium parameters 
BE (b) - The single-medium energy buildup factor 
a,b) - The unscattered energy flux albedo 


2.  Two-medium parameters 


Set oie) The scattered energy buildup factor 


Be (b, 5b.) - The interface energy buildup factor 


ay (b b 


1? >) - The scattered energy flux albedo 
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ieee tb) - The energy flux which scatters in the 
first layer but not the second layer. 
Normalized to an interface flux of 
IMeV/cm*-sec. 
The single-medium parameters are assumed to be available in 
the open literature; however, for convenience, these quant- 
ities are calculated and tabulated in this work for three 
source energies (1, 3, and 6 MeV) and five materials (water, 
aluminum, iron, tin, and lead). The remainder of this 
chapter will be concerned with evaluating the two-medium 
mereameters for two-layered slabs of water, iron, and lead. 
It will then be shown in the next chapter that these quanti- 
ties can be accurately extrapolated to many other two- 


material combinations by interpolating over atomic number. 


4.4.1. Calculation of the Interface Energy Buildup Factor 
and the Scattered Energy Flux Albedo 
The energy buildup factor at the interface between the 
two layers is obtained by applying the energy buildup factor 
definition to the discrete ordinates results for the scat- 
tered energy flux at the interface. Thus, 
el pee bi ) I°>(b,,b 
Ee ieee 2 - 
B-(b,,b,) = 1 + ———— =1+———  . (4.6) 
Bee eliee. 2 U 
This two-medium parameter, as calculated, can't be directly 
extrapolated to other material combination with reasonable 


accuracy; however, accurate extrapolations can be achieved 
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by using the expression for Be (b,,b,) as given by equation 
mo. that is, 

Z Uu S 
B.(b,,b,) - Be(b,) er a (bo) + (B.(b,)-1] a(by>b>) . (4.7) 
As will be shown in the next chapter, the scattered energy 
mux albedo, an (by bo) can be extrapolated to many other 
material combinations. This two-medium parameter is evalu- 


ated for all the two-material combinations of water, iron, 


and lead by rearranging equation (4.7) to yield 


It should be noted that the scattered energy flux albedos, 
as calculated from the above expression, have values ranging 


2 ] 


in order of magnitude from 10 ~ to 10. while the buildup 


Bactors in this equation range in order of magnitude from 


0 to 10°. inicomincdbedtes that third place accuracy in 


10 
the interface and single-medium energy buiidup factors can, 
Bietne average, be expected to provide, at best, first place 
accuracy in the scattered energy flux albedos. At first 
glance, the resulting inaccuracy appears to be a shortcoming, 
however, it must be remembered that the scattered energy 

flux albedos are calculated and tabulated solely for the 


purpose of providing a means to accurately extrapolate inter- 


face energy buildup factors. When equation (4.6) is 
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employed for this purpose, it can be seen that, in reverse, 
first order accuracy in the scattered energy flux-albedo 
provides third place accuracy in the extrapolated interface 


energy buildup factor. 


4.4.2. Calculation of the Energy Flux which Scatters 
mothe First Layer but not in the Second Layer 
In Chapter 3, the energy flux which scatters in the 
first medium but not the second medium, eis ob) was 


normalized and expressed in terms of a truncated Legendre 


expansion of the angular energy flux at the interface be- 





tween the two layers. fhat is, 
5 E S 
O I (b me Ee) 
Mambe+b,) = y Sf) Pee eh) dE 80 (4. 9) 
n - 2 2 cis b.) Q 
g=0 0 |e 


where ib, .b.,£) are the Legendre expansion coefficients, 
Bb, .b.) is the scattered energy flux at the interface, 
(oj) Sas 


mb) = Ete (b) is an exponential integral of 


Q Tie tic 
Mrder jt2, and b = Un (E)X>. Note that b # bo(E,). The 
Legendre expansion coefficients, I°(b,,b,,£) for 9=1,2,...,5, 


gale ae 
are obtained in 20-group structure from the discrete ordin- 
ates calculations while the attenuation coefficient data 
are obtained from the tabulations of Storm and eee! 
These discrete quantities along with the exponential! inte- 


mrals in M, (b) must be expressed as functions of energy in 


order to perform the indicated integration in equation (4.9). 
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A computer program called NINT was developed by the 
author to express the integration variables as functions of 
energy and then numerically integrate equation (4.9). This 
program utilizes a least-squares polynomial-curve-fitting 
routine to approximate and express the attenuation and 


Legendre expansion coefficient data as polynomial functions 


of energy. In formula notation, we have 
J 
Mumeb,.—) = Y g . EN} eM for g=1,29..-.5 
ar Q an 
4 ie see ey 
41-00) 
and 
k 
eee) = x, f oh et, (4.11) 
n=] 


where g and h are the polynomial expansion coefficients, 
and the indices j and k indicate the orders of the polynom- 
jal expansions. (The curve fitting routine was limited to 
M2en Order expansions.) 

The exponential integrals occurring in Mb) are ex- 
pressed as functions of energy by first approximating E, (b) 
and then using a recursion relationship to determine the 
higher order exponential integrals in terms of E,(b). The 


[61 ] 


following Rand approximation for E, (b) is accurate to 


Meeeast One part in a million in the range 1 « b < 110: 


oe >] (4.12) 
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where 
ap = ORS iaeoo Cy = Lede 76/3 3 | 
a 4.5307924 Cy) = 8.6660126 
a5 = 5. 1266902 Cy = Sieooc 7 2 


The higher order exponential integrals are obtained by re- 


peated use of the recursion formula 


I ce . 
a) = —— [e cme) | for n> 1. (4.13) 


With the integration variables expressed as functions 
of energy, equation (4.9) is integrated in NINT with a nu- 
merical integrating subroutine called CADRE. CCeute  oil| 
Adaptive Romberg Extrapolation). The two-medium parameter 
Bb, +b.) is evaluated in this manner and tabulated in 
Appendix C for all of the two-material combinations of water, 
iron, and lead; three source energies (1, 3, and 6 MeV); and 


40 different mean free path combinations. 


4.4.3. Calculation of the Scattered Energy 
Buildup Factor 
The scattered energy buildup factor He (eee) was 
defined in Chapter 3 as a two-medium parameter which gives 


the buildup of the scattered energy flux at the interface as 


* 
CADRE is a CDC 6500 system library program of the 
Computer Sciences Center at Purdue University. 
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jt traverses the second medium. In formula notation, 
; 5 aes) + foto, b,) 
Be(b,,b,) = ————_______+ ,, (4.14) 


Su 
I (b) +b.) 


where i (b, +b.) 1s the energy flux which scatters in both 

layers and seal) has already been evaluated. There is 
no need to evaluate 1°°(b, +b.) because the scattered energy 
buildup factor is simply calculated by rearranging the two- 


medium buildup factor formula to yield 


-b 
Z 
. : [B,(b, +b.) - Be (b,) J e 
[B-(b,.b,)-1] t (b, +b.) 
wtb, .b,), like all of the other two-medium parameters in 


equation (4.15) is evaluated and tabulated for all of the 
two-material combinations of water, iron, and lead; three 
Source energies (1, 3, and 6MeV); and forty different slab 


thicknesses. 


4.5. Determination of the Two-Medium Dose Buildup Factor 
Of all the different buildup factor forms, the dose 


Duildup factor is perhaps one of the more useful to shield 
designers, since many shield systems are designed to limit 
dose rate. In this section, the two-medium dose buildup 
factor will be evaluated by expressing it in terms of the 
previously determined two-medium energy buildup factor. 
This approach has been chosen because it allows the two- 


medium dose buildup factor to be calculated without evaluating 





a new set of two-medium parameters. 
Consistent with the single-medium definition of dose 
buildup, the two-medium dose buildup factor is defined as 


mmemrdtio Of total dose rate to the dose rate resulting 


from the unscattered flux. In formula notation, 
fo ue" (e) 15(b,+b,,£) dE 
meib,+b,) = 1 + SS (4.16) 
D ] 2 fe ee 14 (b +b ) 
O en 0 | ea 


where a, is the energy-absorption coefficient for air, and 


Bb) +b.) is the scattered energy flux spectrum at the 
outer edge of the second slab. Multiplying and dividing the 
numerator in equation (4.16) by the scattered energy flux 


and noting that the unscattered energy flux can be expressed 





as 
Meme hey 8c expl-(b,+b,.)] — oe aa 
1 +5, 3 1755 
Ci SS lae 
yields 
B(b, +b.) =e | + : 
EO va'(e) 15 (b, +b, ,E) 
mE (b.+b,)-1] 1 ne 
Ee | 2 tes [5b +b.) 
jo en O ee 2 
(4.18) 


where B.(b,+b.) is the two-medium energy buildup factor, 
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and Beeb, tb,.£)/1°(b, +b.) is the normalized scattered energy 
flux spectrum which is tabulated in Appendix C. The integral 
in equation (4.18) which represents the average value of 

the quantity we CE) us T(E) is evaluated in terms of the 


20-group energy structures used in the discrete ordinates 


calculations as follows: 


ait 20 aul  e S 

an ; Yen 6E;) 1' (by +bo,E;) 

air iv ) aire ) 7S +b, ) a i 
Fen (E,) 1=] Yon O : (b, 2 


where the average energies E> group widths AE. and group 


* 
averaged energy absorption coefficients of air Pee (G) 


Bapulatced for the 20-group structures of 1, 3, and 6MeV in 


ales 


Tables C.59 - C.61. Substituting equation (4.19) into 


equation (4.18) yields 


7 7 —air, air 
Bib +b) = [B.(b, +b.) a ae [Yon en (4.20) 
Using equations (4.19) and (4.20), the two-medium dose build- 
me, factor, By (b,+b5). and the normalized, average energy- 
absorption coefficient of air, Eee CE.) are calculated 


and then tabulated along with the two-medium energy buildup 


factor in Appendix C. 


* 
The group averaged energy absorption coefficiengs jot 
dir were constructed from the tabulations of Hubbell. 
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[eer neoeNTATION AND DISCUSSION OF RESULTS 


5.1. Description of Tables and Figures 


Almost all of the results obtained in the preceding 
Chapter are presented in tabular form in Appendix C. A Se- 
lected portion of the tabulated information is also reproduced 
in graphical form in Appendix D. It is intended that the 
graphical representations be used primarily to facilitate 
Visualization and analyzation of the extensive collection of 
tables. The purpose of this section is to present the cal oe 
culated results and to aid the reader in wading through the 
mengtchy collection of tables and figures. 

Buildup factors and the parameters in the two-medium 
maeraup factor formulation form the bulk of both the tables 
and figures. For convenience, a summary of the tables which 
defines the tabulated quantities and indicates their location 
and order of arrangement is provided at the beginning of 
Appendix C. The same procedure is employed for the graphical 
representations in Appendix D. In general, the two-medium 
Duildup factors and parameters have been evaluated and tabu- 
lated for all of the two-material combinations of water, iron, 
and lead; three source energies (1, 3, and 6 MeV); and slab 
thicknesses of up to five mean free paths of the first ma- 


terial followed by up to eight mean free paths of the second 
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material. The calculated results for any one of the two- 
medium quantities are grouped first by source energy and 
mien Dy Material configuration. Within this scheme, the 
paired material combinations are arranged by the increasing 
Order of their atomic numbers with the material in the second 
medium taking priority over the material in the first medium. 
For example, application of this arrangement scheme to all 
the paired combinations of water, iron, and lead yields” 

meee liron-water (26 - 6) 

2. Lead-Water (82 - 6) 

3. Water-Iron (6 - 26) 

fe Lead-Iron (82 - 26) 

fee Water-Lead (6 - 82) 

6. Iron-Lead (26 - 82) 
This somewhat unusual arrangement was chosen to emphasize 
mmaeerater utilize the important fact that most of the two- 
medium quantities are similar in both inagnitude and shape 
mereecOntigurations with the same material in the second layer. 

Appendix D presents a small number of tables of Appen- 
dix C in graphical form. These plots are not intended to re- 
place the tables but rather to show how the results change 
With source energy, material configuration, and slab thick- 
Meese being more specific, the reproduction limitations for 


Graph size and grid-line detail are not compatible with the 


“The atomic number of water is taken as six for both 
arrangement and interpolation purposes. 
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medic ions for accurate reading of the graphs. It is there- 
fore anticipated that the reader making serious use of these 
meouits will want to plot his own curves for purposes of 


mmcerpolation. 


5.2. Aids in Applying the Results 


When energy or dose buildup factors for two-layer slabs 
of water, iron, and lead are required, tney can be obtained 
meom the corresponding tables in Appendix © (Tables C.1 - 
meee. Lt will probably be frequently desired to apply 
these results to practical problems which differ in several 
aspects from the calculated results. Source energies may be 
different or the two-layer media may differ in thickness or 
Material from those tabulated. In addition, other buildup 
Tactor forms such as energy absorption may be required. For 
use in such problems, the given results must be interpolated 
Or extrapolated. Some of these inanipulations are quite 
obvious and need no further discussion. For example, the 
Simplest extensions involve obtaining buildup factors for 
different source energies and slab thicknesses. In these 
Cases the tabulated results are plotted and directly inter- 
polated. In contrast, the extension of these results to 
Other material combinations is quite involved because two- 
medium buildup factors can't be directly extrapolated with 
reasonable accuracy. However, as will be shown in Section 
9-3, accurate extrapolations to other materials can be 


achieved by employing the two-medium buildup factor 
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formulation which was derived in Chapter 3. 

When other buildup factor forms are desired they can be 
Obtained in terms of the energy buildup factors by using the 
calculated results for the normalized, scattered energy flux 
at the outer edge of the second slab, Cee) 1 (by tb) 
(Tables C.62 - C.93). For example, consider the energy- 
absorption buildup factor which is defined as the ratio of 
total energy-absorption rate to the energy-absorption rate 


mm~eneresults from the unscattered flux. With this defini- 


meieetne twO-medium energy-absorption buildup factor is 


given as c 
O S 
; Yen te) I (b) +b, >.) dE 
B, (b, +b, ) = ] + en 5h 4.) . Go |) 
Pee. Jie. | ines 7 
en’ oO O 
where Pen is the energy-absorption coefficient for the ab- 


sorbing medium. It should be noted that when the absorbing 
medium iS air, the energy-absorption and dose buildup factors 
are identical. It is therefore not surprising that the pro- 
cedure for evaluating the energy-absorption buildup factor 

Dy expressing it in terms of the energy buildup factor is 
identical to the procedure for obtaining the dose buildup 


factor which was given in Section 4.5. Thus, 


i b,+b 


ae [B. (b,+b,) 7 1 Yen’ Hen \"o 


where the normalized average energy-absorption coefficient 


is given by 
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a 20 


S 
Ven _ ; Men (Ey) 1 (by +b,,E;) AE. (5.3) 
Han ‘Fo! -_ Hen Eo) T° (bye) 


and where the average energies E and the group widths AE . 
meemtabulated for the 20-group structures of 1, 3, and 6 MeV 


in, Tables 359 =- C.6]1. 


oe} Extrapolating the Calculated Results to 


Other Two-Material Combinations 
The two-medium buildup factor formulation was derived 
Mest tOrth in Chapter 3. This formulation will now be 
used to obtain buildup factors for two-material combinations 
other than those of water, iron and lead. For convenience, 


mers rewritten here as 


S SU 
a) = B-(b.) + Be(b,.b,)[B-(b).b,)-1] te (b, +b.) 
- | 2 ‘aes exp Ey 
(5 204) 
where 
u 77.8 
Be (b,.b,) = Be (b,) 2 a (bo) + [B-(b,) 1] v4 (Dy 2b5) (5-5) 


The reason for expecting these formulas to provide accurate 
extrapolations where direct interpolation of the desired two- 
medium buildup factor fails is made clearer by examining the 
two-medium parameters in equation (5.4). In this formula 
there are three, two-medium parameters which must be evaluated 
For meer material combinations before the two-medium buildup 


maetoOr can be determined. These are 
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(1) BE (>, ,b.) - The scattered energy buildup factor, 

(2) i (b, +b.) ~ The energy flux which scatters in the 
first layer but not in the second layer, 
normalized to an interface flux of 1 
Meier seer. 

3 ) Be (b,>b,) - The energy buildup factor at the inter- 


face between the two slabs. 
The first two parameters, BE and are are tabulated in Appen- 
dix C and plotted in Appendix D for all of the two-material 
combinations of water, iron, and lead. The plots show that 
these two-medium parameters are slowly varying functions of 
‘the first layer and smoothly varying functions of the second 
layer. Closer scrutiny of these figures reveals that the 
curves for each of the quantities have the same shape for 
different first layer thicknesses and/or materials. In other 
words, the curves for all the two-layer configurations having 


the same material in the second layer differ only by a cons- 


meant. With this type of functional behavior, it is reasonable 


S 
te 


to other material combinations by interpolating over atomic 


to expect that both B- and ete can be accurately extrapolated 
number. Moreover, it is obvious that these quantities can be 
very accurately interpolated to two-layered slab configura- 
tions with materials in the first layer other than those 
calculated. 

In contrast, the third two-medium parameter, Be (b,,b5) 


1s strongly dependent on the material and thickness of the 
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first layer. Furthermore, this quantity cannot be directly 
extrapolated to other material combinations with reasonable 
accuracy; however, accurate extrapolations can be achieved 

by using equation (5.5). In this formula, the influence of 
the second medium on the interface energy buildup factor is 
represented by the albedos for both unscattered and scattered 
energy fluxes. The unscattered energy flux albedo is nothing 
more than the albedo for a monoenergetic, normally incident 
photon source which has been evaluated in the open litera- 


[5 J 


cure. However, for convenience, this quantity has been 
recalculated and tabulated in Table C.55 for three source 
BEnergies (1, 3, and 6MeV) and five materials (water, alum- 
Baume iron, tin, and lead). During the process of calcu- 
lating the scattered energy flux albedo, it became evident 
that this two-medium parameter could best be characterized 
in terms of the unscattered energy flux albedo for a 

mumpce reflecting medium. Therefore, the ratio, R(b),b,), 
Of scattered to unscattered albedo is tabulated along with 
the scattered energy flux albedo in Tables C.56 - C.58. One 
Mieuie Outstanding features of these results is the rapidity 
with which the albedos assume their maximum value with re- 
Spect to the thickness of the reflecting medium. A re- 
flecting medium of two mean free paths is, for all practical 
Purposes, equivalent to an infinite medium. Hence, these 


albedo results have been tabulated only for reflecting 


medium thicknesses of one and two mean free paths. 
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In formula notation, the ratio of the scattered to the 
unscattered energy flux albedo is given as 
S 
O (b 70 ) 
ay | 4 . (5.6) 
b 
a. | >) 
Miostituting this relationship into equation (5.5) yields 


R( by, 


Be (b) b>) = Be(b,) + a,(by) + [Be (b,) - 1] R(b, bp) a, (bp) 


YY 


c | 
ce 


Mme advantage of extrapolating R instead of a? to other ma- 


terial configurations becomes apparent when one observes that 
this ratio varies only from a value of one to four over the 
entire spectrum of calculated results. This advantage, 

memo led with the fact that the albedo contributions to 

Be (b, »b.) are relatively small, can be expected to yield 
highly accurate extrapolations for the interface energy 


merdup factor. 


cele ine Extrapolation Technique 
Before describing the procedure for extrapolating build- 
up factors with the two-medium buildup factor formulation, 
it is desirable to first define the sets of two-material 
combinations for which the extrapolation technique is applic- 


able. They can most effectively be delineated in terms of 


In Section 4.4.1, it was pointed out that the albedos 
dre approximately two orders of magnitude smaller than the 
BerrespoOnding buildup factors. 
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Z, and Lo» the atomic numbers Of the materials in the first 
and second layer respectively. Using this terminology, the 
range of two-material combinations between water (Z = 6) and 


lead (Z = 82) is given as 


eee 625) © < 25 < 82 


Within the limits of this complete set, eight distinct sub- 
sets are defined and herein considered. These eight ranges 

of material interpolation, as set forth in Table 5.1, have 
Deen cnosen so that the extrapolation technique may be applied 
to and evaluated for most of the two-material combinations 

Of interest in shield design. The extrapolation technique 
Will be described in the remainder of this section and eval- 
Meeead in Section 5.4 by establishing accuracy limits for 

mene of the ranges of interpolation. 

The manipulations and methods employed to extrapolate 
the tabulatcd results can be most easily illustrated by 
Changing the arguments in the two-medium buildup factor form- 
ulation. The mean free path variables in these formulas 
actually represent slab thickness anc material dependence, 
where the material dependence is accounted for with an atten- 
uation coefficient. In this section, slab thicknesses will] 
be held constant and material dependence will be expressed 
in terms of atomic number. Specifically, the independent 
Variables Dy and b. are replaced by Z, and Lo» and the two- 
Seadum buildup factor formulas, equations (5.4) and (5.5), 


are rewritten as 
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e 5.1. Ranges of Material Interpolation. 


Interpolation Range 
imber (in atomic number) 


Pldnomomcerpolation Ranges 


bo < zZ, < 26, 25 = 26 
Zz, = 6, VAS Es Zo < 82 
G < Zz, coe Zo = 82 


—e8e =F ew were ww wy wy www we Wwe we KH wN—— we we Ki wer N' kN’ NX’ NX’ se wX' we ewyew=~ew BF" |\F KF | K| | KF sre eN |g | ZF ee | | — 


Double Interpolation Ranges 


mine. 2. < 82, Gin Zn < 26 


1 


on < 25 20% 26 < Z5 < 82 





and 25 are the atomic numbers of the materials in the 


irst and second layer, respectively. 
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pee (ZZ. )-1 11°" (Ze +7.) 
meetz,) = 8,(Z,) + 5-5, ieee ee 
(5.9) 
and 
Be(Z,.Z,) = Bp (Z,) + u,(Z,) + [B_(Z,)-1] R(ZysZ,)ay (Zp) 


The two-medium buildup factor is extrapolated with equations 


(5.9) and (5.10) by interpolating the calculated results 


j>¥ 


ie and R. More precisely, these 


mor tne quantities Be 
two-medium parameters are linearly interpolated over atomic 
meer With a general set of formulas. For convenience, the 
linear interpolation formulas are presented along with the 
ranges of materia] interpolation in Table 5.2. 

A working knowledge of the above described extrapola- 
tion technique can be achieved by applying the formulas in 
maple 5.2 to a few specific problems. For PAE . consider 
the case of a 3-MeV photon source with four mean free paths 
of tin followed by two mean free paths of water. This two- 
matberia] combination lies within the first interpolation 


mange Of Table 5.2. Therefore, the linear interpolation 


meormula, 


2 Z(Sn)-26 
Met H502) = X(Fe4,H.02) + “Kot £> [x(Pb4,H,02) 
- X(Fe4,H,02)] (5.11) 


(text continued on page 71) 
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Table 5.2. Summary of the Extrapolation Technique. 


I. The Two-Medium Buildup Factor Formula 


Let: 2, = the atomic number of the material in the 
fabieste Layer, 
aa the atomic number of the material in the 
second layer, 


S SU 
Be (Z)5Z5) [B.(Z,,25)-1] ly (2, +Z,) 


exp -Dy 


ee ee ee ee ee ee ee ee ee ee ee ee ee ee) 


meeeecLinear Interpolation Formula for Each Set of Two- 
Material Combinations 


Let X(Z1»Z5) represent any of the following two-medium 
Darameters: 


S 
iP. Be(Z)sZ5) 


SU 
Z Ly (Z,+Z5) 


os U 
S. R(Z,5Z5) = (2, »Zo)/ag (Zo) 


A. Single Interpolation Ranges 


Mezco < Zz, < 82, 


= z(h,0) = 6 


ee 2 


z,-26 
: | : 
2. Za mune) = 26, 6a a= 26 


-6 
eon) = X(26,6) + ee zene?) - X(26,6)] 
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mepie 5.2. Continued. 


] 2 
Z,-6 
(82 Z5) = X(82,6) + SELE reeaeco) = x%( 82,6) | 
4 6 < Zz, < 26, ay = ieee =-26 
-6 
X(Z,,26) = x(6,26) + a (X(82,26) - x(6,26) 
5 a z(H,0) = 6, \a es 82 
Z,-26 
X(6,2Z,) = (6,26) + goog LX(6,82) - es Alo) al 


Mos << 27, « < 26, «222 = z( Pb) = 82 


-6 
;»82) ex (6502) + wie [X(26,82) - X(6,82) ] 


B. Double Interpolation Ranges 


me 26 C< zZ, < Sc. 6 < Zo < 26 
Z5-6 
X(Z1 525) = X(Z, 6) E TELE mrus2-2ove-)x( 82.6) | 


where X(z,,6) is obtained from formula (1) 


mee Cl < CZ, CO «C2, 2b) < 28 fe 


] 


-6 
X(ZsZ5) = X(6,2,) + SEE z bee2o.62) = X(6,82) | 


where X(65,Z5) is obtained from formula (5) 


rr nc rr 








~~ = 
Da 7 






‘ed to interpolate the two medium parameters Be. 


me. Specifically, in formula notation 


2(sn4,H,02) = B2(Fe4,H,02) + 50-26 CBS 


(Fe4, 
E 2 re Pb4,H,02) 


2 ? 


- Be(Fe4,H,02)] (5.12) 


alues for the bracketed quantities are obtained from 


S = 
Be (Pb4,H,02) P22 56 


S = 
ee) Sea eellis 


bstituting these values into equation (5.12) yields 


Aue: <o [2.25 ey 


Be(Sn4,H 02) 


2 


= 2.289 - (5.13) 


reover, the same interpolation formula with 


(Pb4,H,02) = 6.452 x Toe 


(Fed ,H,02) se 


9 pe 6 xX. 1.0 


tained from Table C.45 gives 





i - 
4 - 







0] 2 = 2 


T) 
td 


02) = 5.406 x 10°° + £* [6.45 x 10°° - 5.406 x 10°*] 


Q 
2 


oo. 854 x 10. (5.14) 
ikewise, with the values 
R(Pb4,H,02) = 1.82 and R(Fe4,H,02) = 3.35 
btained from Table C.57, 
R(Sn4,H,02) = 3.35 + £2 [1.82 - 3.35] 
- 2.694 (5.15) 


le interface energy buildup factor for this problem is given 


equation (5.10) as 


eo f(t 
7 \ rem | 14 ’ H 


O02) = B 


2 c ( 


Sn4) + a, (p02) + [B,_(Sn4) 


Aa 


- 1] R(Sn4,H,02 ) 4 


HO (5.16) 


ere the single-medium buildup factor and the unscattered 


Mm albedo obtained from Tables C.54 and C.55 are 


Sn4) = 2.642 and «1°(H,02) = .0249. 


B go Ho 


F ( 


202) pee ee +) .642)(2.694)(.0249) 


| 27 ia Cs ie 
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Now, with all of the two-medium parameters evaluated, the 
two-medium energy buildup factor is calculated with the aid 


Meeeguation (5.9). That is, 





Be(Sn4,H,02)[B,(Sn4,H,02)-1]I°"(Sn4,H,02) 
Re, 
Me ee 229 77) (5.854x10 *) eee 


cee ol Ome 


oc 7 56 


which is the desired result. 
A material extrapolation problem requiring double in- 
terpolation can be illustrated by considering the case of a 


3-MeV photon source with four mean free paths of tin followed 
by two mean free paths of aluminum. This two-material com- 
Dination lies within the seventh interpolation range of 


Table 5.2. Thus, the linear interpolation formula for the 


bee 


é S : . 
two-medium parameters Bes 5 and R 1s given as 


02) + 24Al)-6 [X(Pb4,Fe2)-X(Pb4,H,02) ] 


X(Sn4,A12) = X(Sn4,H ae 


2 
05 V9) 


Where the parameters denoted by X(Sn4,H,02) have been inter- 


Polated in the preceding example. Once more, in formula 


mocation, 


Be (Sn4,A12) = B?(Sn4,H 02)] 


: 02s se O[ BE (Pb4,Fe2)-BS(Pb4.H 


Aicra oosmes 


> 


2 2 


(5.20) 
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| where Be(Sn4,H,02) has been interpolated in the preceding 
example and the values of the bracketed quantities are given 


Mmemiaples C.36 and C.3/7 as 


S 


B°(Pb4,Fe2) = 2.126 and Be 


; (Pb4,H 


292) See 5. 


Substituting these values into formula (5.20) yields 


Le 
20 


B°(Sn4,A12) = 2.289 + 


; RZ ecom- 922564 


2.244 . ore 1) 


Again, the same interpolation formula with eet) = 


Z 


5.854x10 “ obtained from the preceding example and with 


c Ova 6 6052x10 ° ob= 





a 


1°" (Pb4,Fe2) moog and 1°" (Pb4,H 


2 
tained from Tables C.46 and C.45 gives 


2 


5.854x10°*° + ¢ 


i 


154 (Sn4,A12) Pee < = 6.452210 1 


Ma 


S 


2 


= 5a20|1xl0- (522) 


Likewise, with R(Sn4,H,02) = 2.694 obtained from the pre- 


2 
ceding example and with R(Pb4,Fe2) = 1.73 and R(Pb4,H,02 ) 
meee Obtained from Table C.5/7, 
onesie) = 2.694 + a a2 | 
= 2.662 . oeee 3) 


The interface energy buildup factor for this problem is given 


by equation (5.10) as 





q 15 


u 


B°(Sn4,A12) = Be(Sn4) + a, 


: (A12) 


+ [B-(Sn4) - 1] R(Sn4,A12) a (A12) (5.24) 


where the single medium buildup factor and the unscattered 
maperdibedo obtained from Tables C.54 and C.55 are 


Be(Sn4) = 2.642, and a4 (A12) =) 0..023.8 


mnAus , 


me on4),Al2) mee ee 023 +911. 642) (2.66) (0.023) 


c | 


Ze Oo: 5 2 


The two-medium energy buildup factor is now evaluated by 


Meestituting these interpolated values into equation (5.9). 


mat 1S 
Be (Sn4,A12)(B_(Sn4,A12)-1]1°"(Sn4,A12) 
a. 
ee 244 (1.765 5. 88x10 
Hes Ook 1 
a3 698 coer OF) 


Brchn 1S the desired result. 


meee «DISCUSSION Of Accuracy and Comparison of Results 


meorder to use the results of this work it is neces~- 
Sary to have some idea of the expected range of error in 
Doth the calculated and extrapolated results . There are 


Obviously many sources of error in a work of this nature. 
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Some of them are inherent in the computational methods while 
others are due to data inaccuracies. It is convenient to 
massity these sources of error in two categories: 

Meee rrors resulting from incorrectly solving the 
Boltzmann transport equation. 

Meer rors occurring in extrapolating the calculated 
results (especially to other two-material combina- 
moans ). 

In the first category we seek to determine how wel] 

the chosen computational methods solve the transport equa- 
tion by considering the sources of error and comparing the 
meoults. A major goal of this work is to establish defini- 
Meyers accuracy limits for the calculated results. In par- 
mcurdr, it 1S most important to verify the discrete ordin- 
ates results. The importance of this 1S made clear when one 
recalls that most of the parameters in the two-mediuin build- 
up factor formulation were obtained from the discrete ordin- 
ates calculations. Furthermore, the extrapolation technique, 
memset forth in Section 5.3.1, will be evaluated lJater by 
comparing the extrapolated results to corresponding discrete 
Ordinates calculations for a number of two-material combina- 
tions. 

Perhaps the best method of determining the accuracy of 

the discrete ordinates results would be to compare them to 
cOrresponding experimental results. However, for many 


reasons, some obvious, this was not possible. In searching 
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fl 


the literature, the author only found an inadequately small 
number of experimental results for two-layered ilabee oO eae 
On the other hand, if a large number of experimental results 
existed or if they could be easily obtained, there would be 
no need for this work. It goes without question, that an 
undertaking to experimentally verify tne calculated results 
is beyond the scope of this work. 

In the absence of experimental results, it seems rea- 
SOnable to verify the discrete ordinates results by compar- 
ng them to results of other computational methods found in 
mie literature. Here, one finds that the bulk of photon 
attenuation calculations have dealt with single-layered 
Media. However, a number of calculations for multiple lay- 
ered media have been carried out. Some of these of particu- 
lar importance to this work are listed in Table 5.3 along 
With the method of solution. A relevant feature shared by 
mre results of all of these works is that their accuracy 
1S estimated to be at best 10%. A considerable amount of 
computing time and effort has been expended in this work to 
marculate, as accurately as possible, the attenuation of 
Photons in two-layered media. In particular, the discrete 
Ordinates calculations were carried out with the expecta- 
tion of obtaining results which are accurate to within five 
Bercent. It is obvious that the discrete ordinates results 
Bannot be verified to this accuracy by comparing them to 


results which are, at best, accurate to within 10 percent. 
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meee this in mind, the results of the works listed above 
will not be used to establish definitive accuracy limits on 
the discrete ordinates results. Instead, some of them will 
memincluded in the comparisons of Section 5.4.2 with the in- 
mention Of enhancing the credibility of all results in this 
work. 

Without experimental results or calculated results 
which are accurate enough, verification of the discrete or- 
dinates calculations to within 5 percent accuracy is primar- 
ily dependent on comparing these results to the results ob- 
meme rom the stochastic calculations of this work. In 
comparing the results from these two independent methods of 
Solution, it 1s necessary to consider the sources of error 
inherent in each of the methods, paying particular attention 
to tnose sources of error which are conmon to both methods. 
More precisely, errors which are common to both methods of 
Solution must be quantitatively estimated before any mean- 
ingful conclusions can be obtained from the comparison of 


results. 


Meet.  oOurces Of Error in the Calculated Results 

Many possible sources of error are found within both 
Mememods (discrete ordinates and stochastic calculations). 
mous, it may be helpful, before discussing them, to list 
tnese sources of error and to indicate the methods of solu- 


tion in which they are found: 
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Memectatistical errors (stochastic method) 

2. biasing errors (stochastic method) 

3. mesh interval spacing (discrete ordinates) 

4. angular quadrature scheme (discrete ordinates) 

5. group averaging of cross section data (discrete 

ordinates) 
meeetruncation of Legendre expansions (discrete ordin- 
ates) 

Memmmeross section data inaccuracies (both methods) 

eeeeomission of coherent scattering (both methods} 

meeeomission of fliouresence radiation (both methods) 

meee onission of bremsstrahlung (both methods). 
meets section we will first briefly consider the sources 
meer ror which are not common to both methods of solution. 
Mees it 1S Only necessary to show that the resulting inac- 
Curacies are small because the total effect of these errors 
On the attenuation calculations will be accounted for by 
comparing the results of both methods. On the other hand, 
Mes necessary to consider, in detail, those sources of 
error which are common to both methods. Therefore, the bulk 
of this section will be devoted to quantitatively estimating 
mime inaccuracies resulting from common sources of error. 

Stochastic methods, by definition, provide probabilis- 
tic solutions. It is therefore meaningless to present sto- 
chastic results for a quantity of interest without specify- 


Meaethe statistical variance. (Statistical variance is a 





err 


measure of dispersion about the mean value.) In the sto- 
chastic calculations of this work, the variance has been 
estimated in terms of fractional deviation. Specifically, 
a statistical model was developed and used in Section 4.2 to 
express the fractional deviation as a function of the number 
Of histories. With the resulting formulation, equations 
meeiyeand (4.2), the number of histories for each of the 
stochastic calculations was chosen large enough to limit the 
fractional deviation to within 2 percent. It should be noted 
Meeeetne Statistical model, as set forth, does not account 
Moretne effects of biasing. The stochastic code PUGT II 
uses a particle weighting technique to analytically account 
for photoelectric absorption. Pnewsbetisticalevariance 15 
Obviously reduced with this type of biasing. However, the 
resulting decrease in the fractional deviation cannot be 
easily determined. It is the authors belief that the frac- 
tional deviation is less than 1 percent for all of the build- 
Meeeraccor results obtained from PUGT II. In any event, it 
iS known to be less than 2 percent. 

The next four sources of error, numbered three through 
SiX, are pertinent only to the discrete ordinates calcula- 
tions. All of these have been dealt with in Section 4.3 and 


Appendix A. As shown in the discrete ordinates derivation 


* 
, This is the only type of biasing used in the PUGT II 
code. 
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of Appendix A, the definition of a finite-difference phase 
Space mesh and the subsequent approximations involved in 
integrating the differential equations over the mesh consti- 
tute the major sources of error. A great deal of work has 
been done to establish criteria for determining adequate 


[58,65 ,66] The discrete 


angle, energy, and space meshes. 
Ordinates calculations of this work comply with these cri- 
Merida in all cases. 

The error resulting from the group-averaging of cross 
Section data 1s, of course, closely related to the energy 
mesh. It follows that in-group energy dependence can be 
more accurately estimated for finer group structures. The 
multigroup constants, as defined in Appendix A, can be 
exactly determined only when the weighting functions (i.e., 
in-group energy dependence) are exactly known. In practice 
these weighting functions are not precisely known and they 
must be estimated from similar results of other works or 
meeiecalculations with finer group structures. In this work, 
Infinite-medium flux spectra from Reference [1] were used as 
Weighting functions to provide the required multigroup con- 
Stants. However, when the flux spectra, as calculated with 
this group-averaged data, differed significantly from the 
Infinite medium tabulations, the cross sections were re- 
weighted with the calculated flux spectra. 

The remaining approximation which is only relevant to 


the discrete ordinates method is the truncated polynomial 
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mepaneron Of the differential scattering cross section. In 
Meee), the degree of anisotropic scattering dictates the 
order of the Legendre expansion required to achieve accurate 
results. Experience has shown that the most severe aniso- 
tropy, which results from Compton scattered photons with 
energies less than 10MeV, can be adequately resolved with a 
23] 


Pe or Pe 


expansion was used for all of the discrete ordinates calcu- 


expansion. Meiers inamind, a Pe Legendre 
mae1ons in this work. 
The listed sources of error which are common to both 
methods of solution can be categorized as follows: 
1. Omission of photon interaction processes which 
are assumed to have a small effect on the calcu- 
lated results. 
Meeeunaccuracies in the cross section data. 
The subject of photon interactions has been treated at length 


mn the Mierature.-6/ 68.69] 


Many of the more detailed 
treatments indicate that even in the restricted range from 
2O0Kev to 10MeV, the number of processes by which photons 
interact with matter is quite large. It is beyond the scope 
of this work to present an extended description of these 
Various interaction processes. Moreover, it is unnecessary 
Since, as indicated, this area has been excellently covered 
in the literature. However, for the readers convenience, 


Seristing of the various interaction processes which indi- 


cates their relative importance in attenuation calculations 
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Mempresented in fable 5.4. This listing was reproduced 
meome table 5.1 of Reference [27 ]. 

Both the discrete ordinates and stochastic calculations 
meeeinm account for all of the primary interaction processes 
(photoelectric effect, pair production and Compton scatter- 
mma. However, these calculations only account for one of 
the secondary reactions (annihilation radiation). Here, we 
Shall be concerned with quantitatively estimating the effect 
of neglecting the remaining interaction processes which are 
assumed to be of small importance (coherent scattering, 


fluorescence radiation, and Bremsstrahlung). 


Memerent Scattering. In the primary interaction pro- 


cess of Compton scattering, the electrons are considered to 
De unbound. In reality, they are bound to the atom. The 
binding makes little difference in the predictions for 
Compton scattering; however, it does make coherent scatter- 
inG possible. That is the Scattering of photons by the 
Memectron cloud of the atom. In this interaction process, 
the recoil momentum is taken up by the heavy atom. There- 
Tore, the photon experiences almost no change in energy and 
Meeyetittle change in direction. In fact, it has been 
demonstrated that for I1MeV and higher energy photons, deflec- 
tion angles are less than 2° in aluminum and less than 4° 
in Mead. 164 Maproportion to the total absorption coef- 


Ticient, coherent scattering is most important for middle Z 


elements (Z = 25-75) and energies below 500Kev. At higher 





Gamma-ray Interaction Processes 20Kev to 10MeV. 










Of Primary Importance 


Photoelectric effect 
Pair production 
Compton scattering 


Of small importance 


4. Coherent (Rayleigh) electron scattering 
5 Annihilation radiation 

Fluorescence radiation 

Bremsstrahlung 


Negligible 


Thomson scattering from the nucleus 

Beibruck or potential scattering 

Coherent molecular or crystal scattering 
Nuclear interactions 

a. Nuclear photoeffect 

me>- Nuclear scattering 

Radiative corrections to lower order processes 
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energies the cross section falls of rapidly because the 
binding effects are negligible. Further proof of the small 
importance of coherent scattering in shielding calculations 
meeround in Reference {1]. Here, calculations both includ- 
ing and omitting coherent scattering were performed and the 
results differed by only a few percent. 

From the above discussion it seems reasonable to neglect 
coherent scattering. However, it should be neglected in 
the correct manner. The angular distribution of scattered 
photons, as obtained from the Klein-Nishina formula, does 
Mot include coherent scattering. If this interaction pro- 
mss) 1S Chen included in the attenuation coefficient data, 
the resulting effect 1s equivalent to treating coherent scat- 
Bering aS an absorption process which it is definitely not. 
This is the reason that attenuation coefficients without 


coherent scattering have been used throughout this work. 


Fluorescence Radiation. This interaction process re- 
mers to the emission of characteristic x-rays and less ener- 
getic photons by an atom subsequent to the ejection of 
atomic electrons in the photoelectric effect. More precisely, 
X-rays of relatively low energy are given off when free 


electrons fill the electron vacancies created by the photo- 


Memerectric effect. The effects of fluorescence radiation are 


most pronounced for low source energies and high Z materials. 


hOWever, even in these cases, the intensity of this radia- 


meron iS usually Jess than 1% of the intensity of the 











————— 
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radiation producing the original photoelectric effect. In 
fact, according to Reference [70], the fluorescence radia- 
tion contribution to dose rate from a 500 Kev source in lead 
is less than 2%. Hence, it 18 obvious that the omission of 
fluorescence radiation has a negligible effect on the cal- 


mpactved results in this work. 


Bremsstralung. Bremsstrahlung are gamma-rays produced 
when charged particles are decelerated in the atomic electric 
field. The effects of this interaction process are most 
pronounced for high-energy photons (greater than 5MeV) in 
heavy materials such as lead. It has been customary to as- 
Sume that bremsstrahlung are of little importance for source 
energies below 10MeV. However, recent workst? 172! have 
incorporated theoretical models for bremsstrahlung into sto- 
chastic calculations and shown that the effects of this 
interaction process are quite significant for source ener- 
gies of 6 MeV or higher. meron! || has considered in de- 
tail the effects of Bremsstrahlung for photons of 6MeV and 
merey in slab shields of iron and lead. His results clearly 
indicate that inaccuracies resulting from the neglect of 
Dremsstrahlung in the 6MeV calculations of tnis work could 
be as large aS 6% when iron is the material in the second 


layer and 20% when lead is the material in the second layer. 


* 
Recent experimental memisie / have also confirmed 
that the effects of bremsstrahlung are quite significant 
for source energies of 6 MeV -or higher. 





88. 


Fortunately, Dutton has presented his results in such a 
manner that they can be utilized to correct the 6-MeV results 
memes Work. Specifically, a portion of his results are 
presented in terms of bremsstrahlung contributions to the 
total flux spectra (Figures 20, 21, 22, and 25 of Reference 
[71]). With these results as a basis, the following method 
has been devised to approximate the effect of bremsstrahlung 
in the 6-MeV results of this work. 

In what follows, let a prime sian (') denote quantities 
which have been corrected to account for bremsstrahlung and 
let C. denote the correction factor for the energy flux in 
the ith energy group. Then the normalized scattered energy 
flux spectrum at the outer edge of the second layer which 
has been tabulated in Appendix C is corrected to account for 


bremsstrahlung with the formulas 


S > 
(by tb,-€) : 1” (b,+b,,€; ) 


S 1 5 
ib, +b, ) (ole 


aoe tee ee, 208 (5.27) 


where E 1s the average energy of the ith group and the 
Morrection factors C. are obtained from Reference [71]. 


mpecifically, 


Meee When iron is the material in the second layer 


1. C, Peeoeor 0-0 <E. < 1.0 


Ze, C. ic or. | 0 <C. Ores: 


5. C. meer or “il.8 < E 
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Meewnen tead is the material in the second layer 


il. C. =42.0 Tor U0 < E. < eed 


Zz. C. =o for, 1,4 < E a 226 


S. C. Svat Or 20) * E . 


By definition, the bremsstrahlung-corrected two-mediun 
energy buildup factor 1s given as 

| i‘ 15(b,+b,,£) dE 
Mp, tb,) = 1 + —., (5.28) 


b, +b.) 


This quantity is actually evaluated in terms of the 20-group 


mericture Of the 6MeV calculations. That 1s 


! a Rapa. 4) AE. 
a) - 1+ )  — (ees) 
2 14 (b. +b 


1=1 | 2) 


where the average energies ES and the group widths At. are 
Given in Tables C.59 - C.61. Multiplying and dividing 
equation (5.29) by the scattered energy flux leant) and 
noting that 


etperai aD ) 
Be (b, +b.) = | + a ; eseecase 
I (b) +b.) 
wields 
| 20 lee oe.) AE, 
= Be, % ee, ee ne c 
Be{b, +b.) 1 + [B,(b,+b,) 1] (8.3: 


I> (b,+b 


te ] 2) 


where the normalized scattered energy flux spectrum has al- 


ready been corrected to account for bremsstrahlung in 








mauation (5.27). 
In Section 4.5, the two-medium dose buildup factor was 


evaluated in terms of the two-medium energy buildup factor 


as 
—air 
Meb.tb,) = 1 + [B_(b,+b.)-1] oe (5325) 
m | 2 amelie 2 air ; 
Tene) 
where 
mea ir 20 S | 
Ven _ cmc ae 17 (by, +b, .E,) 
—~ .. * SS a a aan 0 [S3c) 
Meo! izt Men'Fo? 1 (by *bo) 


Mrs Gduantity is also corrected to account for bremsstrah- 


lung by again using the bremsstrahlung corrected scattered 


energy flux spectra. Hence, 
“eG aa 
ee) =) Be (bb. )-1) Pee (5.34) 
Heel 2 Roe le 2 Bae) 
where 


a ails S 
air ieee ter) 1 (bo +b, ,£. ) 
ee 5 sir, 4 nea (ORO!) 
air U Gael (bo +b.) 
a (E.) j=) en O eee 


Before estimating the increase in accuracy which re- 
Sults from these corrections, it should be mentioned that 
Dutton’s results are for isotropic sources in single-layered 
Slabs whereas the calculations in this work are for mono- 
directional sources in two-layered slabs. These inconsis- 


tencies can be placed in perspective by referring to the 


20s 
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energy flux spectrum tabulations of Appendix C. In par- 
ticular, it can be seen that the lower energy portions of 
the flux spectra are quite insensitive to the material and 
mapckness in the first layer. Furthermore, the tabulations 
Mereterence [1], indicate that there is not much difference 
meine spectra for isotropic and monodirectional sources. 
It is therefore reasonable to expect the above-described 
method of accounting for bremsstrahlung to provide some im- 
provement in the accuracy of the 6-MeV results. A conserva- 
tive estimate would be that the existing inaccuracies are 
reduced by 50%, thereby limiting the uncertainties resulting 
from the omission of bremsstrahlung in the 6-MeV calculations 
to 3% when iron is the material in the second layer and 10% 
when lead is the material in the second layer. 

ies tinal source of error which is common to both 
Me~tnods Of solution involves inaccuracies in the microscopic 
Meoss section data. fhe status of cross section data is 
mer) Changing, and it is possible that the cross section 


tabulations used in this workl©0>6?] 


may be in error by as 
Much as 5%. This uncertainty obviously results in large 
errors for such quantities as dose rates and energy fluxes. 
For example, at 10 mean free paths, a 5% change in the total 
attenuation coefficient results in a 65% change (e!/%) in 
me unscattered flux. On the other hand, normalized flux 
Spectra and bduildup factors which represent flux ratios 


are known to be relatively insensitive to small changes in 
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mreemicroscopic data. An estimate of this insensitivity 

can be obtained by again referring to Reference [1] where 
calculations both omitting and including coherent scattering 
were performed. For the case of a 4-MeV source in iron, the 
cross section data with and without coherent scattering 
differ by up to 7% for energies above 200 kev. However, 

the corresponding moments method results differed by less 


man 4%. 


5.4.2. Comparison and Verification of the 
Calculated Results 
Both discrete ordinates and stochastic methods have 
been employed in this wprk to calculate, as accurately as 
possible, the attenuation of photons in two-layered slabs. 
meets the purpose of this section to verify the tabulated 
results by comparing dose buildup factors from these two 
independent methods of solution to each other and to cor- 
responding results from the following references: 
meee invariant inbedding calculations of Shimizu, 
Reference [11]. 
meee Monte Carlo calculations of Bowman and Trubey, 
Reference [10]. 
Seeeinpirically titted Monte Carlo calculations of 
Sharp and Carnesale, Reference [12]. 
Before presenting the results of this comparison it should 
be noted that quantities of interest from the above listed 


references are claimed to be accurate to, at best, 10% and 





ke 


that these results represent the more accurate of those 
found in the literature. One of the major goals of this 
work is to establish definitive accuracy limits for all of 
the results. Moreover, the discrete ordinates calculations 
were carried out with the intention of obtaining results 
Which are accurate to within 5%. It is obvious that results 
Which are accurate to, at best, 10% can not be used to ver- 
meye accuracies to within 5%. Therefore, the results from 
the above listed references have been included in the com- 
parisons only for purposes of enhancing the credibility 

Of this work. Without results from other works which are 
accurate enough,definitive accuracy limits will be estab- 
lished by comparing the results from both the discrete 
Ordinates and stochastic calculations of this work. This 
procedure can be expected to yield reliable accuracy limits 
because errors common to both methods of solution have al- 
ready been quantitatively estimated. 

More than 600 buildup factor cOmparisons have been 
carried out to evaluate differences between the ANISN and 
PUGT II results of this work and the corresponding results 
from the above listed references. In general, dose buildup 
factors were compared for two-layered slabs of water, iron, 
and lead; gamma-ray source energies of 1, 3, and 6 MeV; and 
Shield thicknesses of up to five mean free paths of the 
first layer followed by up to eight mean free paths of the 


second layer. The results of these comparisons have been 





——— 








94 


expressed in terms of deviations with respect to the dis- 
crete ordinates dose buildup factors, and summarized in 
Table 5.5 by presenting maximum and average percent differ- 
ences along with the number of comparisons made for each 

of the given configurations. In general, buildup factors 
from this work agreed best with those of Shimizu, Reference 
[11]. As shown in the summary, the maximum differences are 
mmeeless than 10%, the accuracy that Shimizu has placed on 
his results. 

In order to further characterize general trends in 
these comparisons, soine of the more representative results 
meeep!otted in Figures 5.1 through 5.5. In line with the 
Summary, the plots illustrate the consistently good agree- 
ment between the results of Shimizu and this work. More- 
Over, these plots provide some insight to the relatively 
poor agreement between the results of this work and the 


Monte Carlo results of References [10] and [12]. In this 


Case, they indicate that the maximum differences occur at 


the larger mean free path values where statistical approaches 


are known to be less accurate. 


Meee stdbl!ishing definitive accuracy limits it 1s most 


important to measure the differences between the results of 
the discrete ordinates and stochastic calculations of this 


work (ANISN and PUGT II). The comparative summary in Table 


9-5 reveals that dose buildup factors from these independent 


Metnods of solution differed by less than 1.5% for all of 
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Figure 2.) Comparison Of Dose Buildup Factors for 1-MeV 
Photons in Two-Layered Slabs of Water Followed 


by Lead. 
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comparisons. Moreover, it should be noted that over 90% of 
the results differed by less than 1%. These remarkably 
small differences unequivocally substantiate the conclusion 
that sources of error which are not common to both methods 
of sclution are negligible; furthermore, by taking these 
differences into account along with the discussion on common 
Sources of error, the results of this work can be verified 
to the following accuracies. With the exception of the 
6-MeV results, energy flux spectra and buildup factors, as 
tabulated, are conservatively estimated to be accurate to 
Within 5%. The 6-MeV results are known to be less accurate 
| Because both methods of solution do not account for brems- 
Beeeraniung. In fact, it was shown that the inaccuracies re- 
Sulting from the omission of bremsstrahlung could be as 


large as 6% when iron is the material in the second layer 





and 20% when lead is the material in the second layer. Be- 
Cause Of these inaccuracies, a method of approximating the 
Mmerects Of bremsstrahlung was set forth in section 5.4.1, 

) and shown to reduce the existing errors by 50%. When cor- 
rected, the 6-MeV energy flux spectra and buildup factors are 


Meelieved to be accurate to within 


feo When Water is the material in the second layer, 
2. 7% when iron is the material in the second layer, 
and 


3. 12% when lead is the material in the second layer. 








Toe 


mee. 3. Veritication of the Extrapolated Results 

The procedure for extrapolating buildup factors with 
the two-medium buildup factor formulation has been presented 
in Section 5.3 along with the eight ranges of material 
MmcerpOlation. In this section, energy buildup factors for 
a number of representative material combinations are calcu- 
lated with the two-medium buildup factor formulation and 
compared to corresponding results from discrete ordinates 
Calculations. In order to use these comparisons to estab- 
lish reliable accuracy limits, it is necessary to choose 
two-material combinations for which the extrapolation tech- 
meque 1S Known to be least accurate. 

In general, quantities which characterize photon trans- 
port in material media are mainly determined by the varia- 
tion of the attenuation coefficient with energy; moreover, 
Materials of similar atomic number nave attenuation coef- 
ficients which display similar energy dependence. With 
this view in mind, it is reasonable to exvdect that the 
extraplated results will be least accurate for materials 
With atomic numbers which differ significantly from those 
Of water, iron, and lead. For example, consider the fifth 


range of material interpolation, 


2, = Gee2o < Z oc 


2 


Within this range of interpolation, the extrapolation tech- 


maue should be least accurate for the element tin with 
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atomic number 50 because most other materials of interest 
in this range will have atomic numbers closer to either iron 
er lead. 

In line with the preceding discussion, the two material 
ecmbinations which were chosen for the comparative calcula- 
tions have atomic numbers that lie in the middle of the 
interpolation ranges. These representative two-material 
combinations are categorized by applicable interpolation 
range and listed in the comparative summaries, Tables 5.6, 
7, and 5.8. The percent differences, as given in the 
Summaries, represent the maximum deviation between the ener- 
gy buildup factors which were calculated with ANTSN and 
those which were obtained with the two-medium buildup factor 
formulation. In order to characterize general trends, the 
comparative results for two-material combinations in which 
the largest differences occur have been plotted in Figures 
2-6 through 5.9. Both the summaries and tne figures demon- 
Strate that the largest differences occur in the fifth and 
Seventh interpolation ranges where the materials in the 
second layer range in atomic number from 26 to 82. (Even 
in these ranges the differences are at most 5%.) The fact 
that these two ranges of interpolation consistently produce 
the largest differences can be explained by referring to 
Section 5.3 where the functional behavior a; the two-medium 
Parameters has been analyzed. In particular it was shown 
that the two-medium parameters, Be and ae are slowly vary- 


ing functions of the first layer and smoothly varying 
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functions of the second layer; moreover, a detailed examin- 
ation of the plotted results reveals that the curves for 
each of these two-medium parameters have the sae shape for 
different first layer thicknesses and/or materials. Hence, 
meeis reasonable to expect the extrapolated results to be 
least accurate for interpolation ranges which enconpass the 
largest variation in the atomic number of the material of 
the second layer. 

The differences, as set forth in the summaries, clearly 
miagicate that errors resuiting from obtaining buildup fac- 
ors with the two-medium buildup factor formulation are 

1. At most 5% for two-material combinations which lie 

in the fifth and seventh interpolation ranges, 
and 

2. At most 2.5% for two-material combinations which 


iene in all the other interpolation ranges. 





6. CONCLUSIONS AND RECOMMENDATIONS 


eee conclusions 


The theory of photon transmission and buildup in two- 


layered slabs has been investigated and set forth in this 


work. As a result, the two-medium buildup factor nas been 


Sialytically expressed in terms of single medium buildup 


factors and parameters which in turn can be characterized 


ms TuNnctions of atomic number. It should be noted that the 
two-medium parameters in this formulation are not empirical, 
DUt rather they are directly related to the basic physics 

of transport phenomena in two-layered slabs. 

In applying the two-medium builaup factor formulation, 
we assume that the single medium buildup factors and albe- 
dos are available and then seek to determine the two-medium 
Parameters. These quantities have been evaluated by using 
ANISN, a discrete ordinates code, to calculate the attenua- 
tion of photons in two-layered slabs of water, iron, and 
Mead. In the process of tabulating and plotting the results 
it was observed that the two-medium parameters are slowly 
varying functions of the first medium and smoothly varying 
functions of the second medium; moreover, a detailed exam- 
ination of the plotted results (Appendix D) reveals that 


the curves for each of the two-medium parameters have the same 











| shape for different first layer thicknesses and/or materials. 


Because of this functional behavior, these two-medium 
parameters have been accurately extrapolated to a large num- 
ber of two-material combinations by simple linear interpola- 
mon in atomic number. 

Energy and dose buildup factors for more than 15 dif- 
ferent two-material combinations have been calculated with 
the two-medium buildup factor formulation and compared to 
corresponding ANISN results. These comparisons clearly 
Meemeate that the errors resudting from obtaining buildup 
mactors with this formulation are at most 5%. The unique 
feature of the method and results herein described is that 
buildup factors can be accurately hand-calculated using 
fairly simple formulas for most of the two-material combin- 
ations of interest in shield design. In fact, the major 
limitation on the number of materials which can be investi- 


gated is the lack of single medium data. 


6.2. Recommendations 
The methods and results of this work should prove to 
be useful tools of the shield designer, especially in the 
area of determining the optimum two-material combination 
for a given application. Future development and improvement 
should therefore be aimed at extending the application of 
the two-medium buildup factor formulation and thereby in- 


Creasing its usefulness. 










114 


At present the two-medium buildup factor formulation 
has been applied only to normally incident photon sources 
with energies up to © MeV. In order to consider higher 
energies, it 1S necessary to accurately account for the ef- 
fects of bremsstrahlung on the parameters in the two-medium 
buildup factor formulation. This can be accomplished by 
modifying the discrete ordinates code ANISN to include the 
contribution of bremsstrahlung. One possible method of 
including theproduction of bremsstrahlung in a discrete 
ordinates code is to treat this generation process as a 


pseudo-scattering event. In so doing, the angle and energy 


dependent differential cross section for bremsstrahlung can 


be handled in the same manner as the Klein-Nishina formula 


Wfor Compton scattering. 


Extension of the two-medium buildup factor formulation 
to other source configurations such as plane isotropic and 


oblique incident sources is both natural and straightfor- 


mmard. In particular, source configuration differences show 


up as differences in the expressions for the unscattered 


flux terms; therefore, definitions for the two-medium 


Darameters, as well as the methods for calculating them, 


remain the same because these quantities are functions only 


Mr the scattered flux. 


In verifying the extrapolated results, it was shown 


that errors resulting from obtaining buildup factors with 


the two-medium buildup factor formulation are larger for 














material] interpolation ranges which encompass the largest 
'Variation in the atomic number of the material in the sec- 
Marae iayer. Specifically, the largest errors (5%) were con- 
feerstently found to occur in the fifth and seventh interpola- 
tion ranges (see Tables 5.6-5.8) where the materials in the 
second layer, range in atomic number from 26 to 82. The 
error in these interpolation ranges could be reduced by 50% 
by calculating and tabulating two-medium parameters for 
the following material combinations: 

1. Water-Tin 

meee tron-Tin 

3. Lead-Tin 
Meso) doing, errors resulting from obtaining buildup factors 
With the two-medium buildup factor formulation could be 


mumated to 2.5%. 
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APPENDIX A 


A.J. Derivation of the Discrete Ordinates 
) iinans bor t eaemane. 
The general Boltzmann transport equation describing the 


* 
collective particle motion in a unit cell in phase space 


is 


‘where x. is the total macroscopic cross section, o(r, E, 4) 


1s the angular flux in terms of particles per square cm per 





' a per unit energy and per unit solid angle, 

¥V- 8 @ (r, E, &) represents the loss rate due to leakage 
from the phase space cell (dV dE dg), Da o(r,E,%) represents 
the loss rate due to collisions in the phase space cell, 


and §S (r, Be 0) represents the sum of the sources; that is 


Mer, E, ©) = inscattering source + fission source 
+ independent generating source. 
Equation (A.1) is usually obtained by considering the causes 


for particle increases and decreases in the differential 






t 


ithe development given in this appendix will apply to 
both neutron and gamma-ray transport; therefore, the word 
particle will refer to either neutrons or gamma-rays, and 
the symbol £ will be used to signify neutron or gamma-ray 
Cross sections. 








phase space cell dV dE dg In words, the terms of equation 
| 
‘{A.1) represent the following: 


. LEAKAGE + INTERACTIONS = SOURCES. 


‘tn a completely analogous manner, the discrete ordinates 
! 


\formulation of the transport equation can be arrived at by 


considering the events causing an increase or decrease in 


the number of neutrons contained ina finite difference 
| 


| : . 
‘cell. A derivation for one dimensional geometries is given 


‘following the definitions of the coordinate system used and 


the group neutron flux averages. 


A.1.1. Coordinate Systems 
The coordinate systems used in the description of the 
discrete Ordinates equation are in rectangular, spherical, 
M@meecylindrical geometries and are shown in Figure A.1. 
Distances from the origin are given by the variable rs with 
1 representing the spatial mesh interval index for the dis- 
tance r. The cell orthogonal surfaces at distances r; and 


; 


ye are denoted by the terms A. and A. 


47> and the volume 


_ —_ 


enclosed between i and rey, 18 denoted by care the sub- 


script £ denoting the center of the cell. 


Angular distributions are represented with discrete 
rections. The solid angle oP associated with each direc- 
tion is specified by a direction cosine La and a weighting 
function oe which 1S proportional to the surface area sub- 


memaed by the solid angle. 
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more A.|1. Geometrical Functions for One Dimensional 
Geometries. 


Geometry Variable Area Volume 


Plane X. 1 Ke - X. 


i e9 ° z 2 
Cylinder r 2mr. Ta eh) 
phere r amr ors a me 


Energy groups are specified by the group index g with 
the group having the smallest subscript corresponding to the 


largest energy. 


melee. Distribution Function Averages 


The group average angular flux b (7,8) is given as 
oe 
o9(7 8) - | o(r,£,9) dé. (A.2) 
g 
several averages of %, are needed to derive the discrete 
Ordinates transport equation. The mathematical justifica- 
tion for these averages is embodied within the mean value 
veorem of calculus. First, the average of ar over the 
‘volume Veg and the solid angle 2 of the space Cel ES 
2fined as 


Vie ; 
_... 2 jet ray dg. (A.3) 
itd,w g 
> m V. QT 














Then, the average of % over the surface area A. and the 


solid angle 2 is defined as 


] > 


MG.i.m A, w. by (F, 0) dA do (A.4) 
7 om A 7) 

im 
Weal ly, the average of 4 over the cell volume and over the 


Memectional interfaces of the cell is 


erormal definition of the redistribution coefficient a, 
which accounts for the loss of neutrons due to angular re- 
Meeeri bution in curvilinear coordinate systems, wil? be 
given jater. The averages for the group flux es given above 


meerconceptually displayed in Figure A.2. 


A.1.3. Ihe Discrete Ordinates Formulation 
of the Transport Equation 
A flow balance for a finite difference cell wil! be de- 
‘rived in terms of the flux averages given above. The loss 
rate me to leakage from the face of the phase space cell 
ies 


| i | 


iio oO ,1r i,m 





eee 


$g,i+ 1/2, m-1/2 


Pg, i,m Pq i+1/2,m 







g,i+ 1/2 ,m+1/2 





f 


migure A.2. Average Fluxes for Group g. 
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| 
| Equation (A.6) represents the total exit current in the 
direction 2 [ronmmetie face Ae 4] of the cell. The net loss 


by leakage from the phase space cell is thus 


(A A. 6 ). (A.7) 


ead og yitlsm ie Greg iit 


a similar inanner, the net loss rate by leakage from the 


Matcerfaces of the direction cell is 


ays : - Q- : (A.8) 
+3 5Mmts Og yittymts Atlas oo itt m-3 


In the process of a particle traversing a curved (cylinder 


Meeeetie direction of neutron travel is constantly changing. 


j 


which results from this angular redistribution. 


sphere) mesh interval the cosine of the angle assoc 1a tea 
1¢@ a terms in equation (A.8) account for the loss rate 


Mies particle loss rate from collisions in the phase 


“space cell 1s given by the product 


) V b (A.9) 


W * . +) 
in Qed ts 1+4 g,1ti,m 


Where > is the group and space cell averaged macroscopic 


g,it3 
total cross section. Finally, the average source is defined 
Meche Same manner as tne average flux given in equation 
Ss... 2 eo Smeraa tev do . (A.10) 
ett tae Wi g 


Vv. q 
1 : hi) 


The discrete ordinates form of the transport equation 


7 . 
mee ODtained by combining equations (A.7, A.8, A.9, and A.10) 
7 


Wes, 
yield 


. . = + : : 
Wn Yin Ai a ‘Geil sae og sim! Mats. ite Pg, ity meg 


. 2 + W ?. ° V ° 2 
“m-L,i+3 aaron m= £ m “g,itg “itt °g,it+s,m 
= W V o S si e 
fe itt g,i+i.m (oe ey 











* convenience, the centered subscripts are omitted. Then 


fueron (A.11) for each group g becomes 
W uA.) O44) om A. o<) 1 Sines Pm+s 5 Ont Dende 


few» Vd - wl S.~. (Ai) 


The angular redistribution coefficient is associated 
h ray-to-ray transfer in curvilinear coordinates and can 
further described by noting that the net angular redis- 
ibution must be zero. In other words, angular redistribu- 
yn does not create or destroy particles. This condition 

t on a requires that the sum of the terins in equation 
12) over all the discrete directions must equal zero. 
mis, 

MM 

L (Ont Pinta ~ Cm-4 Om- 


] 


OMM+ 4 OMM+ 4 ao! $ = O.. (A.13) 
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In order to meet this condition on the two limit (or end- 
point) values, a and a,, are set equal to zero. Second- 
t MM+2 4 

lly, in the case of divergence free flow (Veto = 0 and 


bea) = d.)> equation (A.12) reduces to Yo = S. Then it 


- A.) (A.14) 


With a 


i 4 for the a's. In the case of plane geometry A. 4, = A. ; 


, = 0, equation (A.14) provides a recursion relation- 
2 


more> equation (A.14) shows that for this case all the 
7 Beare equal to zero. This result follows from the fact 
that there is no angular redistribution in plane geometry. 


does not equal A. 3 however, 


1 : A : 
mv linear coordinates Ae 4] 


follows from equation (A.13) that the m sum of equation 
mia) with the first and last a's equaling zero must vanish. 
ia implies that 


MM 


| i; 6OW = 0 (paeel 53) 


m= | 
f Meation (A.15) represents a constraint which the quadrature 
) tC used to specify the set of discrete directions se must 

Bmeisfy. This is another way of saying that the quadrature 


3 t used must correctly integrrate over the solid angle .. 















A.1.4. Equivalence of the Discrete Transport Equation 
with the Boltzmann Transport Equation 

It is now necessary to show that the discrete ordinates 

orm of the transport equation approaches the analytic form 

miewpoltzmann equation as the finite-difference phase 

ace cell approaches differential size. Dividing equation 


.12) by w V gives 


WAS sy ne eas by) P Om+s?m+s ne ies 
V w V 
+ ¥po = S. (A.16) 


plane geometry as the space and angular mesh intervals 
roach differential size the variables w - dy and V ~ dx. 


mn equation (A.16) becomes 


PG 2S eam 


ch is the form of the Boltzmann transport equation for 
-dimensional plane geometry. 

In spherical geometry as the space and angular mesh 
Mevals approach differential size, w > dy, V > dare aire: 


A > ie The first term in equation (A.16) then be- 


2 2 
; 2 a (ro 9) (A.18) 
ann’ dr re ant 
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relationship for the change in g with respect to di is 


ven by dividing equation (A.14) by w V. The result is 
m+ m-+ _ co: | ie (A.19) 


le limit of this equation as w * du is 


[ews 
& 





ee 
= & e210) 


QJ 


I 
veo. 


tegrating equation (A.20) with respect to uv yields 


2 ee a 
7 al e (A.21) 


meee 1S the constant of integration. Because & vanishes 


the ends of the u interval (a core O) where uw = + 1, 
MM+ 5 


z 
Constant C evaluates to l/r. 
After multiplying equation (A.21) by % and then taking 


Bederivative with respect to u, the resulting expressions 


ad _ | 2 
y = r (] = 4 ) b, 
3 e) 


e left side of equation (A.22) is just the differential 
m of the second term in equation (A.16). Therefore, us- 


Meguation (A.22), it is seen that 


ci mt omy 1 a fer 


es : aS: 
wV ee 6] : 












hus, equation (A.16) becomes 


flee yd = 5. (A.24) 


mise is the analytic (short) form of the Boltzmann equation 
n spherical geometry, where the source term S has yet to 


e evaluated. 


Meleo, EValWation of the Source Term 
The source term 1S in actuality three separate source 


erms as follows: 


S Sei inseattering)/ 


oeitt.m Geet em 


+ S (independent )/ 


g,1ts.m 
aes USS) 000k eae Cae Sy) 
Me analytic form of the inscattering source term is 
| Pees ot) o(r,b 4a’) dE‘ dk J, ene 
E' a | 
(A.26) 


here Ber. E'Es)'+0) Pomeneswscattering cross section for 
hanging the incoming energy and direction £',' to the 
utgoing energy and direction £,{3. The integration is over 
incoming directions and energies. The fundamental 


Ssumption in the discrete ordinates method is that this in- 


scattering integral can be approximated by a quadrature 
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scheme which evaluates the integral as the sum of a discrete 
Meetribution. In one dimensional geometries a Gaussian 
Meddrature scheme can be used to provide the discrete dis- 
tribution, with discrete direction cosines corresponding to 

he Gaussian zeros.. 
The differential scattering cross section can be approxi- 
lated by a truncated Legendre series expansion in the cosine 


Bethe scattering angle one That is 


L 

ta i ' vc a va i 

Bee ES 2) = F Sarto, (FLEE) P (uy), Coop 
=O 


muerte Legendre expansion coefficients, eee ~E)), are 


f 
4, 


iven by 


Myer5E'>E) = peer Sty OQ) PoC.) dy (A.28) 


S L> <0 0 
-] 

Whe scattering angle Uo 1S9vetdted to the initial and final 
Meese and uy by the addition theorem for Legendre 
Miynomials. For one-dimensional plane and spherical 


iometries the relationship is simply 


P (ui) = Fe? Pty) ; (A.29) 
Meeectituting equations (A.27) and (A.29) into expression 
mice) gives 


fete (u) a eee 


Pp a) HE du'/ a iey , oes On) 





Foo 


1e@ multigroup transfer cross section is now defined as 


Eo] ] 
| ceases e)) | Meee ae) Po(u') Ge = cdi). 
g(F.E) = — 
Sigal ; 
eee) Potu') de dips 
EY 


(A.31) 









ubstituting this definition into equation (A.30) and repre- 
enting the integration over E' by a summation of group 


tegrals having the fori 


Etat 
>, (Fou') = Pete sie de CN) 
EY 
elds 
; g | 
Me Pou) F oog(F.€) | opl®u') Pelu') du'/y iain 
h=] -] 


(A.33) 


The generalized Gaussian quadrature formulas, as pre- 
nted in Reference [41], can now be used to approximate the 
tegral in equation (A.33) as the weighted sum of a dis- 


ete distribution. More precisely, 















Bethe direction cosines u. correspond to the roots of 
morenogonal polynomials used to generate the discrete 
ghts ye and the index N represents the order of angular 
idrature (Si); ice wnemmunberwot discrete directions. 
marscrete ordinates form of the inscattering integral is 
Fained by substituting the quadrature formula given by 
gation (A.34) into equation (A.33) and then performing the 
megrations which correspond to the subscripts g,iti, and 


Mie result is 


g N 

kt | h->g c 

2 ae L Ose,its L Wi! P tue) i aT 
h=] m'=] 





APPENDIX B 


Bel. The Relation between the Lattice Model and the 


* 
Integral Transport Equation 


The exact relation between Boltzmann equation and re- 

















urrence relations which define random walks is given by 
Se =. inonul’°! in Appendix 2 of the article’ ®! 
mscussion of the transport equation for particle colli- 
ion density is presented. Here we try to formulate the 
Problem of gamma-ray transport through matter considering 

he probabilistic model on which we have based our sampling. 
It is usual to assume that particle-particle interac- 
tion may be ignored. This is a good assunption in most 
adiation transport problems. It leads to a linear equation 
mamailOws us to discuss one particle at a time in our formn- 
Ulation. The randon walk performed by each photon is de- 
scribed in terms of its collision with atomic electrons and 
macle@i. Between successive collisions it is assumed that 
[the gamma-ray travel in straight line with constant energy. 
Miers further assumed that the expected number of collision 
Meeunit distance is constant. This constant, called linear 
ittenuation coefficient, depends only on the energy of com 
On and the material through which the photon is travelling. 


* 
Reproduced directly from Appendix A of Reference [55]. 





leeks 


yr this analysis we will be concerned only with the steady 













state condition of photon transport, 1.e., time independent 
problem where no transient conditions exist. 
In the model considered, the space is assumed to be 


filled with a cubic lattice and gamma rays are sampled at 


[54 J 


Meee lattice points. iepemoaremdlime together 26 allowed 


directions of motion for gamma rays. We define the location 


of any lattice point by a vector ri ik? the direction of the 


motion of photons by Qe» and its energy by E. The motion 


Mmepnocons can be described by the following quantities. 


aie Pee) eee pected angular density of 


th Collision at 


photons entering their n =a 


per unit volume per unit energy at E, moving 


maedirection, Qn» per unit solid angle. 


Xn i gKe QE) = Expected angular density of 


th 


photons leaving their n eoiision at Poa per 


unit volume per unit energy at E, moving in 


direction, co per unit solid angle. 


t should be pointed out that the energies and directions of 
the x, are in general not the same as for the Vv. 

ie random walk of photon consists of a succession of 
collisions and subsequent traversal of a free path. We de- 


note by Beast: = Pigk? ee = Sao E >» E), the probability 


density that a step changes the state of (Pisgipee Cae) 


mco a unit interval at (755 Mm rE). teas comventent to 


Tactor the above kernel into two functions, one dealing with 





Fog 













changes in the spatial coordinate, the other with changes 


Maeene energy and direction coordinates. 


Me ik mt me +E) = K(r igi pr Ire o E +E) 


en ee 
2 exp [ ag ) MUP pipe tOid Edi] 
q=0 
ee 
6 eee = 2) (A.1) 
gk Magik 


where BG pe tye 2h, OE +E), the Klein-Nishina formula 


meuation (2.2), is the probability at 0 pape of a gamma ray 


‘scattering from ee ) to (2° 4) per unit interval at 
(Q.,£). The exponential function represents simply the 


m 


probability that photon reaches r.., from r,, 


ij at without 
mecollision. d_ is the lattice spacing in the direction 


m 
Ff motion, m. The 6-function represents the fact that pho- 
tons travel in straight line between any two successive 
I11isions, with no change in direction and energy. 

In calculations the sampling is made from the normal- 
Beaeprobability distribution. Let us define the collision 
kernel Cran epi Qe EE) as the expected number of 
yhotons leaving the collision per unit solid angle at es 
Md unit Slew wate 19r a photon entering a collision at 


| eer The relation between the collision kernel 


and the Klein-Nishina formula is 
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me tke tm mE ED (A.2) 


Mora ageypi sk ) 












he first term of equation (A.2) represents the probability 
vat a scattering collision has occurred and the second 
mers the normalized probability distribution for scatter- 


Similarly, we can define the transport kernel 


Be gk mE) as the expected number of the next 


llision per unit volume at ar for a photon leaving a 


meee or a collision at ee te): ee ak: 


»—E) is the normalized transport probability distribution 


d can be found simply by multiplying the second term of 


uation (A.1) by the linear attenuation coefficient, 


rigk?k)- 
re oe 
Me Ke ED = MaKe? Exp[ - ) WIP sr sag 
q=0 
ek 
C Se ee aS 
fae d ,£)d J 6 (— - 2). (A.3) 





14] 


Prom equations (A.1), (A.2), and (A.3) one obtains, 










Mg stk! igke tm nae +E ) = Lu(rasgipeolk ) 
| 





(A.4) 
MCT apoE) 

In a random walk problem we ask in general for the state 

Pi jk? mE) of photon reached after n steps. On the other 

ind the Boltzmann transport equation gives us the total 

melision density, (rs ps2, sE), regardless of the number 

F steps to reach that state. More precisely 


v_( on 


WIrs apoE) a n Pigk em 


K + ©o 


K 
- = Tim ) 


‘It is interesting to note that using Boltzmann theory, we 
ose some information, in that knowledge of Va 7s more in- 
Mnative than the knowledge of V. It should also be 

dinted out that w is a quantity denoting the average be- 
avior of photons in the medium, while from the random walks 
f photons we can obtain information as to the fluctuation 


nd variation from the expected value. Similarly, 


K 

sy laa - = 

X(r 5 p00) ~ 7 ae ) X CP gm E?- 
n=] 
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sare yx and vy, Were previously defined. We note Ke by the 


n 
irce BE ay ED. Considering the definition of y., x,; 
| denoting the lattice volume by v, for photons entering 


4 st , ae - ; 
eir nt] Scattering at position reaps one obtains 


Ue 
mike me? = ) eal cee Or A ane 
7S aa 
25 pr My rE) Av (A.5) 
mee 0, |, 2, 
Similarly, 
KAP gp oo E) ) WyhM gr gr prey sk ) 
m' oie 
Clr Mela sE rE)AQ_, AE ; (A.6) 


By summation over all n, and using the definition of 


W, and yx previously defined, one obtains 


yk? me =) a ) Ue eG yh hE) AV, 


mesimilarly, 


XCr 5 sp 2% E) 


CJ 


VO Km rE’) 
1 F ' 


= 


Crs Meat 7E) AQ dE PS are 















an be eliminated between the equations (A.7) and (A.8), 


obtain 


GK? im? ©? = ) ) Wri gr pr sn sE) 
m 


+02 ,ESE)-T(r. 


me te nt oe AN AV 


* >) etree TAC Deeper, 5, 2 E)Av. (A.9) 


The photon flux, o(r. Pieuemretated to the colli- 


ijko’m? 
yn density by a very simple relation, 

E) = (rs 516) pe ee (A.10) 

The transport equation for flux, 6, can be easily ob- 

ined by substitution of (A.10) into (A.9) and using equa- 
mm A.3), 


ijk mE) i ) ) | ie pe ae se! 
= 


‘ ee a ae 
: ye ake mm? +E )dE AQ 1 AV 
] mm = = = = 
. a ) ee ee a gtk agke mee 


(han) Ty) 
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The integral form of the Boltzmann equation now can be 
meeten by Substituting (A.1) and (A.3) into equation (A.11), 


after summation over the 4-function, as follows 


(5 5,28, 5E) = ) dl ) O(N gp Sd ey 2 MoE ) 
s=0 m' /E! 
° ak om mmm? +E )dE AQ 
S 
+ Sr po Sd 9 8s E) J Exp[- ) ee ds sE)d_ | 
s'=0 
(A.12) 


It is interesting to Arie the similarity of the above 
juation to the familiar form of the integral transport 
Merion <2), except for that in the equation developed 

Dace and direction of motion of photons are discrete quant- 
ities. This is the basic assumption of the model used in 
mers work. In the limit, a lattice distance, d 70 and 

AQ, +0 implying that the number of lattice points and direc- 
10ns become infinite, the two equations should be equiva- 


7 2 
=f rc . 








APPENDIX C 


C.1l. Tables of Results 


The tables in this appendix are arranged as follows: 


Energy and Dose Buildup Factors and 


the Quantity ee ue (ES) for Two- 


Layered Slabs. 


Scattered Energy Buildup Factors, 


23 


p(b,,b,), for Two-Layered Slabs. 


Energy Fluxes, 1 (Basan) for Two- 


Layered Slabs. 


Interface Energy Buildup Factors, 


Be (b,,b for Two-Layered Slabs. 


>) 


Single-Medium Energy Buildup Factors, 
Be (b). 


Unscattered Energy Flux Albedos, 
.U 

& (b 

a4 ( ) 

Scattered Energy Flux Albedos, 


zs 
a,(b,>b5). 


Ci; 


] 


3.0 


42 


seo 


04 


eee. 


eis) 


Gr .c 


C4) 


C50 


Ceo 


Oe Xs! 











ass Energy-Absorption Coefficients 
of Air in the Energy Group Structures 


of the Calculated Results. 


meeecered Energy Flux Spectra, 
A (b,+b,,£.)/1 (b, +b.) for Two 


Layered Slabs. 
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moe C.33. 


for IMeV Photons 


Scattered Energy Buildup Factors, Be 


Iron-Water and Lead-Water. 


E 


(by. 


in Two-Layered Slabs of 


bo). 








3 (0) slag) Iron-Water 1MeV 
ah 
b. 2 3 4 5 6 7 
meoon 2.158 2.165 2.169 2.154 2.176 2.198 
2 Bee 3323 3.348 3.366 3.357 3.397 3.435 
3 4.670 4.678 4.713 4.748 4.720 4.750 4.780 
4 mieten) 43° 6.218 6.243 6.223 6.260 6.290 
5 mer 707.82 7.878 7.937 7.898 7.920 7.940 
6 9.472 9.550 9.666 9.859 9.831 9.871 9.905 
7 bieo ieee) 1),.77 12.00 11:96 12.00 12.03 
8 Meee ieee 14.20 14.35 14.25 14.29 14.33 
Sets ce) Lead-Water 1MeV 
a 
b. 2 3 4 5 6 7 
meioeee 20S) 1.962 1.937 1.897 1.865 1.833 
2 Seoeec 074. 2.998 2.963 2.899 2.867 2.838 
3 Mceeeareil) 49279 4.229 4.110 4.060 4.010 
4 6.364 5.889 5.724 5.659 5.607 5.540 5.493 
5 7.956 7.647 7.423 7.333 7.148 7.060 7.004 
6 9.993 9.392 9.309 9.287 9.201 9.080 8.980 
7 Pee i.66° 11-53 11.40 11.29 11.20 
8 eee ooeia ig 14.12 13.99 13.88 13.79 


a at sy SS >= SS a gap fifi s-e 








maple C.34. 





peaetenea Energy Buildup Factors, B2(b,,b 


(ine 
for 1MeV Photons in Two-Layered Slabs of 


a)» 


Water-IJron and Lead-Iron. 


(b, ,b Water-Iron 1MeV 


as 2) 


err eects oe Il .964 1.963 1.950 1.939 
meow ooo ee. 639 2.840 2.848 2.814 2.796 
Sees 72 02/95 3./94 «83.796 3.770 3.746 
4.795 4.820 4.839 4.835 4.848 4.826 4.798 
peepee oe O00 56 y9e) 6 ©656-99306CU6 005 5.980 5.954 
emer ace eit) fo. 2/2 06/«41 301 ~667.270 «37.240 
633 beoommconol0 8610/7 6./20 8./00 8.6/0 
Vere 0e-15 


TGEab3 Os 


2D) Lead-Iron 


Pesto? ss laoo5 |F65| 1.844 1.835 1.825 
Peeeeces io 2-/07 2.683 2.674 2.664 2.655 
eee ooo oly S262 3.584 3.574 3.564 


seem oeeloemoeoll 9.771 5.752 5.740 5.735 
eee oOe 7 OO) 7.167 7.127 7.115 7.100 


BS es 


On 29 


8.726 


LeOi29 


ey 516 a0 
POR22 


180. 





ros 








Maple C.35. Scattered Energy Buildup Factors, ye (le Ne 
for 1MeV Photons in Two-Layered Slabs of 
Water-Lead and Iron-Lead. 
S 
Be(b, .b,) Water-Lead 1MeV 
] Z 5 4 5 6 7 


eee i 0S) 1.40) 1.394 1.385 1.365 1.346 


meee l.oet 1.885" 1.867 1.852 1.840 
eee ic. Woe 2.093 2.083 2.068 2.058 
eects) 25299) 2.282 2.268 2.248 2.240 
Peoomece= 475 2.49606 2.456 2.444 2.420 2.410 
merece O27.) 2.609 2.590 2.570 2.560 
Meme OG ec. /060 2./39 2./18 2./00 2.690 


»b,) Iron-Lead 1MeV 


mercer 1.4048 1.395 1.374 1.367 1.360 
eooeeeeoet) 1.656 1.647 1.643 1.628 1.605 
mevlcmmeocome|.097 1.886 1.869 1.852 1.840 


meee esos e.cGd 2.275 2.264 2.240 2.230 
moO omece4 ian 2.965 2.45) 2.439 2.418 2.408 
pee oe ols) 2.599) 2.584 2.566 2.556 
eee eee See. fer 2.705 2.680 2.6/0 














mapie C.36. 
bs i 
2 2 
3 2 
4 3 
5 3 
6 4 
7 5 
8 5 
D 
= 
2 2 
3 2 
4 3 
5 3 
6 4 
7 4 
8 5 


ihe 2 


Scattered Energy Buildup Factors, Jeep: 


for 3MeV Photons 


in Two-Layered Slabs of 


Iron-Water and Lead-Water. 


S 

Be(b,»b,) Iron-Water 3MeV 

] 2 3 4 5 6 
aes ies cl Il.7/33 $|1.746 1.748 #1 
auc eee) eel Ss  2.32606-62.328 «= =O? 
molec wo4o = c.ooe 2.908 2.909 2.910 2 
meee oo 4oo | 6 S2500 «63.509 3.510 3 
Piveeeeuos 4.076 4.5103 4.113 4.120 4 
peeeeoce b> 0 646677 «648. /06 064712 C4 
Wem eee e264 5.290 5.299 5.310 5 
Beemer ooo 59 6 6«U856 FOB Yl6hlCUCOGUl BIB C6 .908) OCOD. 

S 

Br (b,.b,) Lead-Water 3MeV 

2 3 4 5 6 

Preece ooy  1o66/ 1.674 1.679 #1. 
meee Clie cee 39) | 24256 2.2/3 %2.2//7 2 
ali~eecen Gece. O3506CU2 O63 6CUZ«. B70 = 02 B74 2 
Pape o moet ob)6U6O8459 3.478 3.482 3. 
See ceo o027,)6 6h 060:hCU 081) OU4 086 4. 
ame ema n bb) 4.650 4.689 4.698 4. 
meow oely oo 23G 5.281 5.30] 5. 
ewer omen 765° 5.659 5.894 5.920 5. 








maple C.37. Scattered Energy Buildup Factors, 52 (Dae) 
for 3MeV Photons in Two-Layered Slabs of 


Water-Iron and Lead-Iron. 





Be (b, 5b.) Water-Iron 3MeV 
a 

b. 2 3 4 5 6 7 
fe Mey eeeieco  ise24 1.644 1.638 1.639 1.640 
1g Ecce zeiet 2014) 2.173 2.167 2.162 2.158 
LE Pememe 68S 2.708 2.753 2.728 2.720 2.710 
Ly 4 ie oe?ss acs 3.345 3.316 3.302 3.280 
(= 3.758 3.875 3.896 3.965 3.923 3.908 3.883 
BR 6 Wesati472 4.523 4.590 4.539 4.521 4.495 
| 7 Moco oO) 50158 5.240 5.184 5.160 5§.135 
1® PcewSn752 55825 5.935 5.875 5.840 5.810 
| Babb, Lead-Iron 3MeV 
| 

b, 2 3 4 5 6 7 
7, Mea o6en 1582 1.593 1.600 1.604 1.608 

2 MOON enO7o 2009 2.126 2.143 2.145 2.146 
ik Poe G6) 24706 2.712 2.714 2.715 

4 ees 275 3.291 3.308 3.313 3.317 

5 37667 3.820 3.879 3.899 3.916 3.925 3.930 

6 f2s494.413 4.485 4.520 4.555 4.568 4.578 

7 emer 5 08) 5.119 5.141 5.153 5.163 


8 ee eomoo oe 7 9255 5./52 5.771 5.784 5.794 





PO oer 
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| meee ©C.38. Scattered Energy Buildup Factors, Bate, PS) 
for 3MeV Photons in Two-Layered Slabs of 
Water-Lead and Iron-Lead. 

) Beton) Water-Lead 3MeV 

| b 

| b., 1 2 3 4 5 6 7 

/ ] eee imo cces les2o 1.332 1.330 1.328 

1% 2 1.575 1.587 1.601 1.605 1.614 1.610 1.606 
3 meso oemioco sia g02. 1.91) %1.910 1.905 1.900 
4 emcee 2-255 2.216 2.223 2.216 2.209 
5 Mees 26556 2.042 2.547 2.533 2.522 
6 eee osc 4Gool 2.662 2.875 2.860 2.848 
y Sev oeeseou ssn 17/ 3.184 3.202 3.192 3.182 

| 8 eetoomeortoy oe l0 S.530 3.550 3.540 3.530 

| Beato D)) Iron-Lead 3Me V 

| 3 

b., ] 2 3 4 5 6 is 

] iweeeeeowo ns i2ocs Tass) 1.336 1.335 1.334 

| 4a Peewee le o0y  le613 1.621 31.615 1.610 

3 Porte osaeie 900 I2916 %1.917 1.913 1.909 

| 4 elomecetioleeecaclo 2.221 2.227 2.2\/7 2.210 

| 5 eee eA noe 2.543 °° 2.551 2.541 2.532 

| 6 eromee eo oem ce 2G 00 2c. 659 2.865 2.85/ 2.852 
/ eens acemooly// oe lot 3.198 3.188 3.178 
8 See oor 16 jemoz6 63.545 34515 3.525 





Table €.39. Scattered Energy Buildup Factors, Sle Ald th 
for 6MeV Photons in Two-Layered Slabs of 
Iron-Water and Lead-Water. 
32 (4 abe) Iron-Water 6MeV 
a 
b, ] 2 3 J » 6 7 
] meet coee le o86 8 12520 1.530 1.533 1.535 
2 Pmovomelsoo5 1.09) 1.908 1.920 1.912 1.910 
3 femesomece2os8 2-202 2.300 2.310 2.308 2.305 
4 metoomee oe). 2.667 5 2,686 2.692 2.694 2.696 
D ere O07 oneo/) 3056  3.06/ 3.073 3.080 
6 eos cd5 so. 575 3-398 3.414 3.420 3.425 
7 Pepe oo 3-094 3.720 3.734 3.739 3.744 
8 portiees- 245 3.979 4.003 4.020 4.025 4.030 
Slay eley Lead-Water 6MeV 
* 
ba, ] 2 3 4 5 6 7 
ii meme tog o75) 1.4696 1.517 1.525 1.531 
le 1.720 1.809 1.854 1.884 1.911 1.909 1.904 
3 Peete eee ces) 2.262 2.291 2.291 2.292 
4 euler of ee 2.500) 2.624 27.650 2.655 2.660 
5 eee ecco ee O96 68.991) 62.997 )~ = 3.005 
6 fee memoria 2073-2498 3.2/8 3.284 3.290 
/ pee ore es 4767 °53.520) 3.553 3.558 3.565 











Table C.40. 


b 

| 
2 
3 2 
4 2 
5 2. 
6 3 
7 3 
g 4 

b 

> 
2 
3 2 
4 2 
5 2 
6 3 
7 - 
8 4 


SPeauccreacs energy Buildup Factors, Be (b 


for 6MeV Photons 


Water-Iron and Lead-Iron. 


Be (b,»b 
1 
vs oe 
mes (| 
mie. 2 
554 2 
988 3 
459 3 
no47 4, 
442 4. 
S 
Be(b,.b 
1 
Po26. 
m666 1. 
moe 2. 
Rileie 2. 
815. 2. 
Rot Se. 
669 3 
ee 


>) 


404 
60 
aU ite: 
ler’ 


Water-Iron 


-942 
409 


LP 


in Two-Layered Slabs of 


1.808 
Fae 


A263 6 


443 
O29 


P2018) 


6MeV 


2.7248 


2.680 


4.548 


| RAL rt 


4.581 


186 








maple C.41. 


ae) 


CO N Do Oo Re | PO 


Scattered Energy Buildup Factors, B2(b,,b 


zl | 
for 6MeV Photons in Two-Layered Slabs of 


>) 


Water-Lead and Iron-Lead. 


Be (b ey) Water-Lead 6MeV 

] 2 3 4 5 6 fe 
.179 .204 med acc 22 | ee Tae 
383 .409 23 yeas 0 Ae miles ey 
.615 . 660 .673 .679 .680 _ oer Dbesys. 
. 883 926 ~942 .950 952 956 958 
159 mec yc 39 pea? S256 Keo4 Pes) 
.468 Boo boo foo .569 Fee ee: 576 
OS .8/78 898 .906 .910 .914 918 
eia2 MOLC nc 43 -254 258 BAO / 266 

S 
Be (b »b,) Iron-Lead 6MeV 

] Z S 4 5 6 7 
ly 4 eelOe -221 wooo eo os EOS mea 
376 416 pss naa 448 ~45) 155 
oy 656 BOWS .689 .696 OS 710 
4 Soy 926 94/7 .96 1 .96/7 ~979 990 
lea Bee lec 238 254 L268 MCT 5 5 gene’ 
~-449 522 oro .569 Poo | Poo 604 
774 855 883 .903 917 .928 e057 
.- 109 .194 acco ese o2'63 me P/aoks 


[icon 
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Energy Fluxes, (yeh). for 1MeV Photons in 


Two-Layered Slabs of Iron-Water and Lead-Water. 


SU 
i (b) +b.) 


] 


.392(-1) 
.706(-2) 
.093(-2) 
.362(-3) 
.063(-3) 
.419(-4) 
.116(-4) 
.681(-5) 


SU 
Ly (b,+b,) 


.597(-1) 
.441(-2) 
.345(-2) 
.202(-3) 
.340(-3) 
.343(-4) 
.424(-4) 
mle? (5 ) 


1MeV 


1MeV 


Iron-Water 
2 3 4 
554(-1) ome = 1) 638(-1) 
144(-2) .301(-2) 368(-2) 
218(-2) peo =2 ) 274(-2) 
HESTON es 8) .841(-3) B75 > 3 ) 
174(-3) .205(-3) 211(-3) 
762(-4) .849(-4) 858(-4) 
224(-4) .248(-4) 248(-4) 
024(-5) .096(-5) 086(-5) 
Lead-Water 
2 3 4 
.919(-1) .098(-1) .202(-1) 
.498(-2) .100(-2) 27a = 2) 
.681(-2) Oud (ac } .000(-2) 
.278(-3) .917(-3) .310(-3) 
.688(-3) .897(-3) MUI OA = 3 ) 
.479(-4) .165(-4) .585(-4) 
.798(-4) .026(-4) eee) 
.960(-5) .724(-5) .178(-5) 


.682(-3 
.148(-3 
.992(-4 


.301(-1) 
.814(-2) 
124(-2) 


eee eee “eee” 


.300(-4) 
.634(-5) 








Energy Fluxes, ia 
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b, tb,), for 1MeV Photons in 


Two-Layered Slabs of Water-Iron and Lead-lIron. 


Su 
is (b.+b 


] 
] 


.374(-1) 
-691(-2) 
.088(-2) 
.352(-3) 
m59(-3) 
.410(-4) 
.113(-4) 
mo? 1(-5) 


im (b.+b 
nN 


] 


.738(-1) 
.944(-2) 
.497(-2) 
.678(-3) 
.492(-3) 
medio (-4 ) 
.585(-4) 
.249(-5) 


Water-lIron 


2 3 
.445(-1) .480(-1) 
.912(-2) 3.964(-2) 
mlss(=c) 1.158(-2) 
.536(-3) 3.536(-3) 
.113(-3) 1.109(-3) 
.570(-4) 3.542(-4) 
.161(-4) 1.149(-4) 
"819(=5) 3.768(-5) 

Lead-Iron 

2 3 
038(-1) lees | =|) 
937(-2) 6.421(-2) 
816(-2) 1.976(-2) 
704(-3) 6.227(-3) 
825(-3) 1.997(-3) 
923(-4) 6.487(-4) 
944(-4) lsiz (= 4)) 
442(-5) 7.068(-5) 


1MeV 


4 


.472(-1) 
O36 = 2) 
.148(-2) 
.485(-3) 
.089(-3) 
.468(-4) 
.122(-4) 
.671(-5) 


1MeV 


.264(-1) 
.746(-2) 
.085(-2) 
mao 3) 
mie (=3) 
.879(-4) 
.262(-4) 
5 03(-5.) 


.464(-1) 
.892(-2) 
.129(-2) 
.417(-3) 
.064(-3) 
.381(-4) 
.091(-4) 
5640-58) 


325 — le 
.972(-2) 
eG ey 
.839(-3) 
.198(-3) 
.157(-4) 
.359(-4) 
.814(-5) 








Table C.44. 


| 





SU 
ie | 


Energy Fluxes, ie 


Ee 


b, tbo), for 1MeV Photons in 


Two-Layered Slabs of Water-Lead and Iron-Lead. 


] 
| 


m0) (-1 ) 
.470(-2) 
.016(-2) 
.109(-3) 
.792(-4) 
.147(-4) 
.027(-4) 
.396(-5) 


f(b +b 


] 


.375(-1) 
me 76{-2) 
m7 7 (-2) 
.299(-3) 
.039(-3) 
.342(-4) 
.091(-4) 
.609(-5) 


Oo 


8. 
ae 


oye 


Water-Lead 1 MeV 
2 3 4 

/ HONDLS) .282(-1) .248(-1) 1 
.465(-2) .389(-2) .279(-2) 3 
.009(-2) 622(-3)} .457(-3) 9 
Or On=3:) .982(-3) .860(-3) 2 
.655(-4) .329(-4) .918(-14) 
.094(-4) .982(-4) .843(-4) 
.007(-4) .689(-5) ei o:) 
.324(-5) ol 5) POsy(=5) 2 

ITron-Lead 1MeV 

2 3 4 

.450(-1) .465(-1) .463(-1) 1 
.863(-2) .886(-2) pac Oae) 3 
.126(-2) ma 2 ) mies =2) 1 
.434(-3) .426(-3) mo76(-3) 3 
.078(-3) OV 2(=3) .053(-3) 1 
456 { -4 ) .429(-4) esemt=4) 3 
.126(-4) .114(-4) .091(-4) 1] 
.714(-5) .674(-5) Poe/(-5) 3 


eon) la) 
.182(-2) 
Si (4G) 
2754(=3) 


563(-4) 
724(-4) 
819(-5) 


.895(-5) 


.453(-1) 
.815(-2) 
mo (= 2) 
e311 = 38) 
030-39 
.281(-4) 
.063(-4) 
.490(-5) 








lwepre C.45. 


Energy Fluxes, 1 ( 


Lack) 


b, +b.) for 3MeV Photons in 


Two-Layered Slabs of Iron-Water and Lead-Water. 


Su 
(b,+b,) 


] 


.791(-1) 
.024(-2) 
.495(-2) 
.586(-3) 
.437(-3) 
.582(-4) 
.478(-4) 
.837(-5) 


3MeV 


Iron-Water 
2 5 4 
See e937(-1) 1.947¢-1) 
Bee oe sog(-2) 5.406(-2) 
.577(-2) oe G2) BoeG(-2) 
.818(-3) .847(-3) .820(-3) 
p55 (3) .509(-3) .497(-3) 
.779(-4) .803(-4) .758(-4) 
.537(-4) 5536(-4) .518(-4) 
.017(-5) 5.004(-5) 4.939(-5) 
Lead-Water 
2 3 4 
.176(-1) .224(-1) ey es ) 
mia (=2) .366(-2) 452(-2) 
.837(-2) .888(-2) 908(-2) 
moe = 3 ) .760(-3) 807(-3) 
.756(-3) .795(-3) 806(-3) 
.610(-4) .693(-4) 709(-4) 
.795(-4) .827(-4) 835(-4) 
.858(-4) 5.952(-4) gear -4)) 


.944(-1) 
.399(-2) 
.581(-2) 
.785(-3) 
.482(-3) 
.681(-4) 
.505(-4) 
.870(-5) 








— 
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Maple C.46. Energy Fluxes, eee) for 3MeV Photons in 
Two-Layered Slabs of Water-Iron and Lead-Iron. 
i (b. +b Water-Iron 3MeV 
] 2 S 4 2 
marat-1) 2.015(-1) 2.046(-1) mg] ) .061(-1) 
Mev /(-2) 5.868(-2) 5.948(-2) .969(-2) .965(-2) 
meo7(-2) 1.790(-2) 1.810(-2) .811(-2) .807(-2) 
meioe-3) 5.609(-3) 5.655(-3) .650(-3) $627 (=3) 
.717(-3) 1.790(-3) 1.802(-3) .797(-3) .788(-3) 
.609(-4) 5.819(-4) 5.844(-4) .809(-4) 1767(=45 
.833(-4) 1.902(-4) 1.908(-4) .898(-4) .893(-4) 
mest=-5) 6.300(-5) 6.311(-5) .271(-5) e216 (= 5p) 
® (b, +b Lead-Iron 3MeV 
———— 
] 2 3 4 5 

Met -1) 2.353(-1) 2.419(-1) 4521 = 1) .470(-1) 
mite) 6.932(-2) 7.127(-2) moe c= ) .264(-2) 
966(-2) 2.119(-2) 2.175(-2) .198(-2) p20 ay 
179(-3) 6.637(-3) 6.799(-3) moo7 4-3) .869(-3) 
Py ie=s) 22119{(-3) 2.164(-3) ge = 3) .179(-3) 
419(-4) 6.860(-4) 6.993(-4) .029(-4) .013(-4) 
PO8(-4) 2.247(-4) 2.295(-4) .307(-4) .294(-4) 
990(-5) 7.440(-5) 7.559(-5) Boot 5 ) nove | = 5m) 
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EMergy Fluxes, Tg ig for 3MeV Photons in 


Two-Layered Slabs of Water-Lead and Iron-Lead. 


SU 
i. (b,+b 


] 


>) 


.339(-1) 
.248(-2) 
.318(-2) 
N55 74 = 3) 
.499(-3) 
.360(-4) 
.823(-4) 
.608(-5) 


Water-Lead 
2 3 
.283(-1) 2.305(-1) 
.099(-2) 7.161(-2) 
Bye) 2290 (=2) 
427(-3) 7.465(-3) 
460(-3) 2.468(-3) 
231(-4) 8.250(-4) 
782(-4) 2.784(-4) 
A72(-5) 9.472(-5) 
Tron-Lead 


.364(-1) 
.309(-2) 
Bose (-=2) 
.584(-3) 
.503(-3) 
moo (-4 ) 
.818(-4) 
.579(-5) 


3MeV 


.311(-1) 
ple Ga } 
.290(-2) 
Boe eats.) 
.460(-3) 
.216(-4) 
.771(-4) 
.417(-5) 


3MeV 


.369(-1) 
.310(-2) 
.326(-2) 
B52 (= 3.) 
.489(-3) 
.298(-4) 
.794(-4) 
.486(-5) 


.311(-1) 
159s) 
.285(-2) 
s422(=39 
.448(-3) 
.168(-4) 
.752(-4) 


.348(-5) 


.368(-1) 
e283 =ee 
.315(-2) 
.504(-3) 
.469(-3) 
.224(-4) 
.769(-4) 
Po Si =i5e) 
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Table C.48. Energy Fluxes, I°"(b)+b,), for 6MeV Photons in 
Two-Layered Slabs of Iron-Water and Lead-Water. 
| Ppa +5, ) Iron-Water 6MeV 
lf 2 3 4 5 
| 
Fa Meo) 9531-1) 2.009(-1) 2.029(-1) 2.033(-1) 
| 2 5.089(-2) 5.529(-2) 5.649(-2) 5.669(-2) 5.647(-2) 
3 Me509(-2) 1.630(-2) 1.657(-2) 1.655(-2) 1.642(-2) 
4 4.596(-3) 4.945(-3) 5.005(-3) 4.981(-3) 4.925(-3) 
5 Meic0(-3) «1.533(-3) 1.547(-3) 1.535(-3) 1.514(-3) 
6 4.523(-4) 4.836(-4) 4.868(-4) 4.821(-4) 4.745(-4) 
7 1.451(-4) 1.548(-4) 1.555(-4) 1.537(-4) 1.511(-4) 
8 4.713(-5) 5.017(-5) 5.031(-5) 4.966(-5) 4.873(-5) 
> +b.) Lead-Water — 6MeV 
| b 
bs 1 2 3 4 5 
1 Mere 2016(-1) 2.086{-1) 2.110(-1) 2.112(-1) 
2 5.076(-2) 5.656(-2) 5.798(-2) 5.805(-2) 5.759(-2) 
3 Mero?) 15659(-2) 1.687{-2) 1.677(-2) 1.651(-2) 
4 Ms5s(-3)) 5.010(-3) 5.066(-3) 5.007(-3) 4.901(-3) 
5 Mets) 1.548(-3) 1.559(-3) 1.533(-3) 1.495(-3) 
6 4.463(-4) 4.873(-4) 4.887(-4) 4.790(-4) 4.654(-4) 
7 MeO ey 1557(-4) 1.557(-4) 1.521(-4) 1.474(-4) 
8 Meee) 5.037(-5) 5.024(-5) 4.898(-5) 4.735(-5) 
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526 (=a 
119(-3) 
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695(-4) 


.888(-4) 


671(-5) 


09 =m 
585(-2) 
349(-2) 


.419(-3) 


380(-3) 
736(-4) 
542(-4) 


Energy Fluxes, I°"(b,+b,), for 6MeV Photons in 
Two-Layered Slabs of Water-Iron and Lead-Iron. 
eb. +5.) Water-Iron 6MeV 

Gg 2 3 4 
1 Mee32(-1) 2.427(-1) .498(-1) Posii=)) 2. 
2 meo74(-2) 7.612(-2) Wocme= 2.) 02) °°) | 274) a a 
3 me2e(-2) 2.422(-2) .489(-2) Sy 2 (9 nas 
1 m7 4(-3) 7.808(-3) .016(-3) .081(-3) 8. 
q mos 4{-3) 2.548(-3) Mo (3 ) moeecl=3) 2. 
6 7.740(-4) 8.401(-4) .608(-4) .664(-4) 8. 
/ mony 6(-4) 2.795(-4) .862(-4) .879(-4) 2 
8 mmpo5(-5) 9.375(-5) eo G25) .643(-5) 9. 
Pb, +b.) Lead-Iron 6MeV 

>) 

b, 1 2 3 4 
] meodO(-1) 2.334(-1) .441(-1) semi -1) 2. 
2 mee/i(-2) 7.168(-2) .460(-2) poead(-2) 7. 
3 1.982(-2) 2.255(-2) .334(-2) .354(-2) 2. 
4 meco7(-3) 7.211(-3) masiht= 3.) .465(-3) 7 
5 mpou7y2(-3) 2.338(-3) .401(-3) maOs=3) 2. 
6 6.820(-4) 7.669(-4) .850(-4) eoos(-4) 7. 
/ 2.266(-4) 2.541(-4) .593(-4) mocu(-4) 2. 
8 7.587(-5) 8.487(-5) neat =5 | Meck -5) 8. 


435(-5) 
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364 eee 2 ieee 3 he 204 eee oo ioe O4 
1.698 eee 2 lee S 2 ees 7 7 1.41] ee lal 
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Table C.55. Unscattered Energy Flux Albedos, ay (b). 


Material Energy U U U 
MeV a (1) a (2) a, (~) 
] 0.0890 .0994 1002 
Water 3 020238 0.0249 oe.02 50 
6 One si2 Oiie3 3 0133 
] 0.0840 (Osis 0920 
luminum Gr0iZ22 .0229 0230 
6 0.0164 .0166 0166 
| OF 06S0 0697 0701 
ron 5 0.0183 0.0184 0.0184 
6 Oh Oa1ONS 0208 0208 
] 0.0280 .0290 0290 
n 3 2010s .0104 0104 
Oe 0) rai 0188 0188 
] Oe eo 0110 
ad 5) Oe Bais 0055 0055 
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Table C.58. Scattered Energy Flux Albedos, a (b).b,), for 


v) 
6MeV Photons in Two-Layered Slabs. 
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a(b,.b,) O41 o 4 a4 (by >bo) O41 04 
b = | b,=2 b,=1 b,=2 b,=1] b= b,=] Dee 
] me 5/ 0.0464 3.46 3.49 b20295 0.0305 2.24 lp tues, 
Z mea 77 0.0482 3.61 3.62 Deez 4 0.0323 2.45 2.49 
3 meee 706 60.0475 3.56 SS Haves 0.0336. 2.50 BOS 
4 meoos 0.0468 3.51 B22 Osa 0 03405 22.53 2.56 
5 meys «60.0477 3.58 3.59 Os0s29 070345 2.57 2.59 
7 meooe 0.0507 3.80 tye OO so Oteie 0 3658 ec ../2 2.74 
Water-Iron Lead-Iron 
] Meoa50 0.0468 2.20 2 ue O20 ei a0) 034 4ae 1.59 6s 
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3 meugo8 0.0476 2.28 2.29 O20 3360-03847. 264 1263 
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fF mm433 0.0441 2.11 Zee e034 SO 0S52 2 |. 6/7 1.69 
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3 meo2z0) 0.0201 1.30 30 Due leo Oe 1 oO Tele Vale 
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Table C.59. The Mass Energy-Absorption Coefficient of Air 
in the Energy Group Structure of the Calculated 


Results for a |]MeV Source. 


Group E. (MeV) AE. (MeV) air 


(7) cm? / gm 
ie 
0.995 0.01 0.0278 
2 0.945 0.09 0.0281 
3 0.850 0.10 0.0286 
4 0.750 0.10 0.0291 
5 0.650 0.10 0.0294 
6 0.550 0.10 0.02955 
7 0.450 0.10 0.02955 
8 0.350 0.10 0.0292 
9 0.250 0.10 0.0278 
10 0.175 0.05 0.0259 
1 15125 0.05 0.0241 
12 0.095 0.01 a 0236 
13 0.085 0.01 0.0241 
14 0.075 0.01 0.0258 
45 0.065 0.01 0.02895 
| 16 0.055 0.01 0.03555 
| 17 0.045 0.01 0.05375 
18 0.035 0.01 0.10745 
| 19 0.025 0.01 0.3295 
| 
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Table €.60. The Mass Energy-Absorption Coefficient of Air 
in the Energy Group Structure of the Calculated 


Results for a 3MeV Source. 


Group E. (MeV) AE. (MeV) u alee 


ay emo /gm 
p 

2.995 0.01 0205 
2 2.895 0.19 0208 
3 2.700 0.20 02135 
4 2.500 0.20 02195 
5 2.300 0.20 02254 
6 2.100 0.20 0231 
7 1.900 0.20 0238 
8 1.700 0.20 0246 
9 1.500 0.20 .0254 
10 1.300 0.20 02635 
1] 1.100 0.20 0273 
12 0.900 0.20 02835 
13 0.700 0.20 0292 
14 0.500 0.20 02955 
15 0.300 0.20 0285 
16 0.150 0.10 0250 
17 0.090 0.02 02385 
18 0.070 0.02 02735 
19 0.050 0.02 04465 
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Table C.61. The Mass Energy-Absorption Coefficient of Air 
in the Energy Group Structure of the Calculated 


Results for a 6MeV Source. 


Group E. (MeV) AE. (MeV) u air 


(—=*) em*/ gm 
0 

1 5.995 0.01 0.0164 

2 5.695 0.39 0.0166 

3 5.400 0.40 20170 

i 4 5.000 0.40 0.0174 
im 5 4.600 0.40 0.0179 
6 4.200 0.40 0.0184 

7 3.800 0.40 0.0190 

8 3.400 0.40 ORO 19.5 

9 Be000 0.40 0.0206 

10 2.600 0.40 OWvaly, 
1] 2.200 0.40 0.0228 
2 1.800 0.40 0.0242 
13 1.400 0.40 0.0259 
14 1.000 0.40 0.0278 
15 0.700 0.20 0.0292 
16 0.500 0.20 0.0296 
47 0.300 0.20 0.0285 
18 OLS 0 o. 0 02 0250 
19 0.080 0.04 O02 56 


20 0.040 0.04 O13 7 
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Table C.84. Scattered Energy Flux Spectra, I°(b,+b,,£,)/ 
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APPENDIX D 


D.1. Graphical Display of Results 


The figures in this appendix are arranged as follows: 


Energy Buildup Factors, Be (b,+b,) for Deal se Dates 
Two-Layered Slabs. 

Scattered Energy Buildup Factors, De 19-0) eee 
eet, +b.) for Two-Layered Slabs 

Energy Fluxes, Ui CO ae for Two- ie Caan hole 


meeayered Slabs 
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Figure D.1. The Energy Buildup Factor, B,(b,+b,) vs. b, and 
b, for 1MeV Photons in Slabs of Iron Followed 
by Water. 





WATER 
LEAD - WATER, b, = 
LEAD - WATER, b, = 


LEAD - WATER, b, 


E = |] MeV 


0 


B. (b, + b,) 





| 2 3 4 9 6 — 68 S, lo We le 


b,+b,, TOTAL THICKNESS (mfp) 


megure D.2. The Energy Buildup Factor, Be (b, +b.) vs. b) and 
b. for 1MeV Photons in Slabs of Lead Followed 


by Water. 
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Figure D.5. The Energy Buildup Factor, Be (b,+b,) vs. b, and 
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Figure D.12. The Energy Buildup Factor, B.(b,+b,) vs. Bb, and 
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Figure D.16. The Energy Buildup Factor, Be (by +b.) VS. b, and 
b, for 6MeV Photons in Slabs of Lead Followed 
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Figure D.20. The Scattered Energy Buildup Factor, Bl baibes 
vs. b, and bs for 1MeV Photons in Two-Layered 
Slabs of Water-Iron and Lead-Iron. 
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Figure D.21. The Scattered Energy Buildup Factor, B  (oepey 
VS. b; and bs for 1MeV Photons in Two-Lavered 
Slabs of Water-Lead and Iron-Lead. 
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megure D.22. The Scattered Energy Buildup Factor, Baba bey 
So by and b., for 3MeV Photons in Two-Layered 
Slabs of Iron-Water and Lead-Water. 
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Figure D.24. The Scattered Energy Buildup Factor, BE (Da abey 
VS. Dy and D. for 3MeV Photons in Two-Layered 
Slabs of Water-Lead and Iron-Lead. 
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Figure D.25. The Scattered Energy Buildup Factor, Dalal.) 
VS. b, and b» for 6MeV Photons in Two-Layered 


Slabs of Iron-Water and Lead-Water. 
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Figure D.26. The Scattered Energy Buildup Factor, Beles aly! 


VS. b, and D, for 6MeV Photon in Two-Layered 
Slabs of Water-Iron and Lead-Iron. 
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mpoure D.2/7. The Scattered Energy Buildup Factor, eS) side) 
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Figure D.30. The Energy Flux, 1°" (b,+b5) vs. b, and by for 
1MeV Photons in Two-Layered Slabs of Water- 


Lead and Iron-Lead. 
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Figure D.32. The Energy Flux, eae) vs. b, and by for 


3MeV Photons in Two-Layered Slabs of Water- 
Iron and Lead-Iron. 
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Figure D.33. The Energy Flux, (Deep VS. b, and b, for 
3MeV Photons in Two-Layered Slabs of Water- 
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migure 0.34. The Energy Flux, 1°" (by +bo) vs. b) and b, for 
6MeV Photon in Two-Layered Slabs of lron- 
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Figure D.35. The Energy Flux, I°°(b,+b,) vs. b, and b, for 
6MeV Photons in Two-Layered Slabs of Water- 
Iron and Lead-Iron. 
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Figure D.36. The Energy Flux, Crete vs. by and by for 
6MeV Photons in Two-Layered Slabs of Water- 


lead and Iron-Lead. 
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